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1.0 Configuration Determination 

To determine the best configuration from all candidate configurations, it was necessary 
first to calculate minimum system weights and performance To optimize the design, 
necesary to vaiy configuration-specific variables such as total system weight, thrust t 
SSSSLSjto durations, total thrust available, and mass fraction for the system. 
Optimizing each of these variables at the same time is technically unfeasible and not 
necessarily mathematically possible. However, discrete sets of data 
which will eliminate many candidate configurations. From the most promising 
remaining designs, a final configuration can be selected. 

To ease this optimization process an interactive spreadsheet 

thrust to weight ratios, number of main engines, vehicle diy weight, andjnitial (p e 
separation) burn time. From this information, it calculated maximum weights 
propellant mass fractions available for each vehicle and the final post-separation bum 

time for the orbiter. 

After selecting the number of engines and thrust to weight ratios for both the booster and the 
orbiter the length of the initial bum was varied to maximize fuel available to the orbiter for 
the second phasebum. This maximum weight was then used to determine tank sizings 
which dictated total volume necessary within the booster and orbiter. 

Included are the three meet important designs considered, oh .which eta*'■W™^*** 
the design criteria set forth in a Marshall Space Flight Center study of the Shuttle II, the 
configuration used in the initial proposal; and the final configuration presented an 
analyzed in this report Note that these three systems represent important iterative steps in 

our final design. 

A listing by cell of the formulas used to generate the aforementioned data is included for 
reference. 
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ORBITER DAT* - 


.Thrust/Weight 

1.2000 

Initial Boost Engines 

o 

o 

Total Thrust at Takeoff 

1, 454,100.00 lbf 

Available Takeoff Weight 

1,211,750.00 lbm 

Dry Weight 

180,000.00 lbm 

Payload 

60,000.00 lbm 

RCS Fuel (MMH/N04) 

20,000.00 lbm 

Total Orbiter Weight 

260, 000.00 lbm 

Maximum Fuel Weight 

951,750.00 lbm 

Fuel Available Initially 

951,750.00 lbm 

Initial Boost Engines 

4.00 

Initial Fuel Consumed 

505,080.80 lbm 

Primary Fuel Available 

446,669.20 lbm 

Auxilliary Fuel 

496,952.57 lbm 

Total Fuel Available 

94 3, 621.. 7 7 lbm 

Final Boost Engines 

4.00 

Approximate Final Thrust 

1,740,000.00 lbf 

Approximate Final T/W 

1.4456 

BOOSTER DATA - 


Thrust/Weight 

1.4000 

Initial Boost Engines 

5.00 

Total Thrust at Takeoff 

1, 817,625.00 lbf 

Available Takeoff Weight 

1,298,303.57 lbm 

Dry Weight 

160,000.00 lbm 

RCS Fuel (MMH/N04) 

10,000.00 lbm 

, Total Booster Weight 

170,000.00 lbm 

Fuel Weight 

1, 128,303.57 lbm 

Fuel Available Initially 

1,128,303.57 lbm 

Initial Boost Engines 

5.00 

Initial Boost Fuel 

631,351.00 lbm 

Reservoir Fuel Pumped 

496,952.57 lbm 

OVERALL DATA ------- 


Thrust/Weight 

1.3034 

Initial Boost Engines 

9.00 

Total Thrust at Takeoff 

3,271,725.00 lbf 

Available Takeoff Weight 

2,510,053.57 lbm 

Dry Weight 

340,000.00 lbm 

Payload 

60,000.00 lbm 

RCS Fuel (MMH/N04) 

30,000.00 lbm 

Total Dry Weight 

430,000.00 lbm 

Fuel Weight 

2,080,053.57 lbm 

Fuel Available Initially 

2,080,053.57 lbm 

Initial Boost Engines 

v> 

o 

o 

Initial Boost Fuel 

1,136,431.80 lbm 

Fuel Available Finally 

943, 621.77 lbm 


Initial Boost Burn Time 
Final Boost Burn Time 
Total Burn Time 

Tsl 

Tvac 

Wdot 

SG - Hydrogen 
SG - Oxygen 
Density of Water 

Max Fuel - Hydrogen 
Max Fuel - Oxygen 
Volume - Hydrogen 
Volume - Oxygen 
Total Volume 
Eq. Hydrogen Tanks 
Eq. Oxygen Tanks 
Minimum Length Necessary 


Max Fuel - Hydrogen 
Max Fuel - Oxygen 
Volume - Hydrogen 
Volume - Oxygen 
Total Volume 
Eq. Hydrogen Tanks 
Eq. Oxygen Tanks 
Minimum Length Necessary 


133.00 sec 

248.48 sec 

381.48 sec 

363 , 525.00 lbf 
435 , 000.00 lbf 

949.40 Ibm/sec 
0.0700 ft~3 
1.1490 ft *3 
62.43 lbm/f t A 3 

135,964.29 lbm 
815, 785.71 lbm 
31,112 .40 ft*3 
11,372 . 68 ft^3 
42 , 485.07 ft~3 
88.0301 ft 
21.4203 ft 
109.4504 ft 


161,186.22 lbm 
967,117.35 lbm 
36,883.88 ft"3 
13,482.36 ft^3 
50,366.23 ft * 3 
69.4704 ft 
25.3939 ft 
94.8643 ft 
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ORBITER DATA - 


Thrust/Weight 

1.2000 

Initial Boost Engines 

5.00 

Total Thrust at Takeoff 

1, 817,625.00 lbf 

Available Takeoff Weight 

1,514,687.50 lbm 

Dry Weight 

180 , 000.00 lbm 

Payload 

60,000.00 lbm 

RCS Fuel (MMH/N04) 

20,000.00 lbm 

Total Orbiter Weight 

260 , 000.00 lbm 

Maximum Fuel Weight 

1,254,687.50 lbm 

Fuel Available Initially 

1,254,687.50 lbm 

Initial Boost Engines 

5.00 

Initial Fuel Consumed 

564,893.00 lbm 

Primary Fuel Available 

68 9, 794.50 lbm 

Auxilliary Fuel 

563,410.57 lbm 

Total Fuel Available 

1,253,205.07 lbm 

Final Boost Engines 

4.00 

Approximate Final Thrust 

1,740,000.00 lbf 

Approximate Final T/W 

1.14 99 

BOOSTER DATA - 


Thrust/Weight 

1.4000 

Initial Boost Engines 

5.00 

Total Thrust at Takeoff 

1,817,625.00 lbf 

Available Takeoff Weight 

1,298,303.57 lbm 

Dry Weight 

160,000.00 lbm 

RCS Fuel (MMH/N04) 

10,000.00 lbm 

Total Booster Weight 

170,000.00 lbm 

Fuel Weight 

1,128,303.57 lbm 

Fuel Available Initially 

1,128,303.57 lbm 

Initial Boost Engines 

5.00 

Initial Boost Fuel 

564,893.00 lbm 

Reservoir Fuel Pumped 

563,410.57 lbm 

OVERALL DATA - 


Thrust/Weight 

1.2923 

Initial Boost Engines 

10.00 

Total Thrust at Takeoff 

3,635,250.00 lbf 

Available Takeoff Weight 

2,812,991.07 lbm 

Dry Weight 

340,000.00 lbm 

Payload 

60,000.00 lbm 

RCS Fuel (MMH/N04) 

30 , 000.00 lbm 

Total Dry Weight 

430 , 000.00 lbm 

Fuel Weight 

2,382,991.07 lbm 

Fuel Available Initially 

2,382,991.07 lbm 

Initial Boost Engines 

10.00 

Initial Boost Fuel 

1,129,786.00 lbm 

Fuel Available Finally 

1,253,205.07 lbm 


Initial Boost Burn Time 
Final Boost Burn Time 
Total Burn Time 

Tsl 

Tvac 

Wdot 

SG - Hydrogen 
SG - Oxygen 
Density of Water 

Max Fuel - Hydrogen 
Max Fuel - Oxygen 
Volume - Hydrogen 
Volume - Oxygen 
Total Volume 
Eq. Hydrogen Tanks 
Eq* Oxygen Tanks 
Minimum Length Necessary 


Max Fuel - Hydrogen 
Max Fuel - Oxygen 
Volume - Hydrogen 
Volume - Oxygen 
Total Volume 
Eq. Hydrogen Tanks 
Eq. Oxygen Tanks 
Minimum Length Necessary 


119.00 sec 
330.00 sec 
449.00 sec 

363 , 525.00 lbf 
435 , 000.00 lbf 

949.40 ibm/sec 
0.0700 ft A 3 
1.1490 ft A 3 
62.43 lbm/ft~3 

179,241.07 lbm 
1,075,446.43 lbm 
41,015.32 ft * 3 
14,992.55 ft * 3 
56,007.87 ft^3 
116.0496 ft 
28.2383 ft 
144.2879 ft 


161, 186.22 lbm 
967,117.35 lbm 
36,883.88 ft *3 
13,482.36 ft*3 
50,366.23 ft*3 
69.4704 ft 
25.3939 ft 
94.8643 ft 


3 






ORBITER DATA - 


Thrust/Weight 

1.2000 

Initial Boost Engines 

4.00 

Total Thrust at Takeoff 

1, 454,100.00 lbf 

Available Takeoff Weight 

1,211,750.00 lbm 

Dry Weight 

180,000.00 lbm 

Payload 

60,000.00 lbm 

RCS Fuel (MMH/N04) 

20,000.00 lbm 

Total Orbiter Weight 

260,000.00 lbm 

Maximum Fuel Weight 

951,750.00 lbm 

Fuel Available Initially 

951,750.00 lbm 

Initial Boost Engines 

4 . 00 

Initial Fuel Consumed 

562.044.80 lbm 

Primary Fuel Available 

389,705.20 lbm 

Auxilliary Fuel 

544,897.09 lbm 

Total Fuel Available 

934,602.29 lbm 

Final Boost Engines 

4.00 

Approximate Final Thrust 

1,740,000.00 lbf 

Approximate Final T/W 

1.4566 

BOOSTER DATA -- 


Thrust/Weight 

1.4000 

Initial Boost Engines 

6.00 

Total Thrust at Takeoff 

2, 181,150.00 lbf 

Available Takeoff Weight 

1,557,964.29 lbm 

Dry Weight 

160,000.00 lbm 

RCS Fuel (MMH/N04) 

10,000.00 lbm 

Total Booster Weight 

170,000.00 lbm 

Fuel Weight 

1,387,964.29 lbm 

Fuel Available Initially 

1,387,964.29 lbm 

Initial Boost Engines 

6.00 

Initial Boost Fuel 

843,067.20 lbm 

Reservoir Fuel Pumped 

544,897.09 lbm 

OVERALL DATA - 


Thrust/Weight 

1.3125 

Initial Boost Engines 

10.00 

Total Thrust at Takeoff 

3,635,250.00 lbf 

Available Takeoff Weight 

2,769,714.29 lbm 

Dry Weight 

340,000.00 lbm 

Payload 

60,000.00 lbm 

RCS Fuel (MMH/N04) 

30,000.00 lbm 

Total Dry Weight 

430,000.00 lbm 

Fuel Weight 

2,339,714.29 lbm 

Fuel Available Initially 

2,339,714.29 lbm 

Initial Boost Engines 

10.00 

Initial Boost Fuel 

1,405,112.00 lbm 

Fuel Available Finally 

93 4,602.29 lbm 


Initial Boost Burn Time 
Final Boost Burn Time 
Total Burn Time 

Tsl 

Tvac 

Wdot 

SG - Hydrogen 
SG - Oxygen 
Density of Water 

Max Fuel - Hydrogen 
Max Fuel - Oxygen 
Volume - Hydrogen 
Volume - Oxygen 
Total Volume 
Eq. Hydrogen Tanks 
Eq. Oxygen Tanks 
Minimum Length Necessary 


Max Fuel - Hydrogen 
Max Fuel - Oxygen 
Volume - Hydrogen 
Volume - Oxygen 
Total Volume 
Eq. Hydrogen Tanks 
Eq. Oxygen Tanks 
Minimum Length Necessary 


148.00 sec 

246.10 sec 

394.10 sec 

363 , 525.00 Ibf 
435 , 000.00 lbf 

949.40 lbm/sec 
0.0700 ft *3 
1.1490 ft ^3 
62.43 lbm/ft~3 

13 5, 964.29 lbm 
815,785.71 lbm 
31,112.40 ft^3 
11,372 . 68 ft ^3 
42,485.07 ft*3 
88.0301 ft 
21.4203 ft 
109.4504 ft 


198,280.61 lbm 
1,189,683.67 lbm 

45.372.10 ft A 3 

16.585.10 ft~3 
61, 957 .20 ft~3 

85.4579 ft 
31.2379 ft 
116.6958 ft 
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'' ORBITER DATA - 

Thrust/Weight 
Initial Boost Engines 
Total Thrust at Takeoff 
Available Takeoff Weight 
Dry Weight 
Payload 

RCS Fuel (MMH/N04) 

Total Orbiter Weight 
Maximum Fuel Weight 

Fuel Available Initially 
Initial Boost Engines 
Initial Fuel Consumed 
Primary Fuel Available 
Auxilliary Fuel 
Total Fuel Available 
Final Boost Engines 
Approximate Final Thrust 
Approximate Final T/W 

BOOSTER DATA - 

Thrust/Weight 
Initial Boost Engines 
Total Thrust at Takeoff 
Available Takeoff Weight 
Dry Weight 
RCS Fuel (MMH/N04) 

Total Booster Weight 
^ Fuel Weight 

Fuel Available Initially 
Initial Boost Engines 
Initial Boost Fuel 
Reservoir Fuel Pumped 

OVERALL DATA - 

Thrust/Weight 
Initial Boost Engines 
Total Thrust at Takeoff 
Available Takeoff Weight 
Dry Weight 
Payload 

RCS Fuel (MMH/N04) 

Total Dry Weight 
Fuel Weight 

Fuel Available Initially 
Initial Boost Engines 
Initial Boost Fuel 
Fuel Available Finally 


1.2 

5 


-B3*E6 

lbf 

-B4/B2 

lbm 

180000 

lbm 

60000 

lbm 

20000 

lbm 

-B6+B7+B8 

lbm 

-B5-B9 

lbm 

-B10 

lbm 

-B3 


—B13*18*X2 

lbm 

-B12-B14 

lbm 

— B3 4 

lbm 

-B15+B16 

lbm 

4 


-B18*E7 

lbf 

—B19/(B9+B17) 


1.4 


5 


—B23 *Z6 

lbf 

-B24/B22 

lbm 

160000 

lbm 

10000 

lbm 

-B2 6+B27 

lbm 

—B25-B28 

lbm 

— B2 9 

lbm 

-B23 


— B3 2 *18*12 

lbm 

-B31-B33 

lbm 

-(B4+B24)/(B5+B25) 


-B3+B23 


— B4+B2 4 

lbf 

-B38/B36 

lbm 

-B6+B26 

lbm 

1 

0> 

lbm 

-B8+B27 

lbm 

•B40+B41+B42 

lbm 

•B10+B29 

lbm 

-B12+B31 

lbm 

-B37 


■ B47 *E2 *E8 

lbm 

-B46-B48 

lbm 
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Initial Boost Burn Time 
Final Boost Burn Time 
Total Burn Time 

Tsl 

Tvac 

Wdot 

SG - Hydrogen 
SG - Oxygen 
Density of Water 

Max Fuel - Hydrogen 
Max Fuel - Oxygen 
Volume - Hydrogen 
Volume - Oxygen 
Total Volume 
Eq. Hydrogen Tanks 
Eq. Oxygen Tanks 
Minimum Length Necessary 


119 

sec 

-B17/(B18*E8) 

sec 

-E2+E3 

sec 

363525 

lbf 

435000 

lbf 

949.4 

lbm/sec 

0.07 

ft A 3 

1.149 

ft A 3 

62.43 

lbm/ft*3 

-B12/7 

lbm 

-B12*6/7 

Ibm 

-E13/(E9*X11) 

ft*3 

-E14/ (XI0*E11) 

ft A 3 

-E15+E16 

ft *3 

-X15/ (PI() /4*15 A 2) /2 

ft 

-E16/ (PI()/ 4*2 6 A 2) 

ft 

-E10+X19 

ft 


Max Fuel - Hydrogen 
Max Fuel - Oxygen 
Volume - Hydrogen 
Volume - Oxygen 
Total Volume 
Eq. Hydrogen Tanks 
Eq. Oxygen Tanks 
Minimum Length Necessary 


B29/7 

lbm 

B2 9*6/7 

lbm 

127/(19**11) 

ft * 3 

128/(110*111) 

ft A 3 

129+130 

ft *3 

129 / (910/4*26*2) 

ft 

130/(PI()/4*2fi A 2) 

ft 

132+133 

ft 
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2.0 STME Performance Calculation Summary 

A model of the total thrust available from the Space Transportation Main Etipnes used on 
rtw. Shuttle II was necessary to evaluate its launch profile. Thrust is a function of many 
S^ indST^Kansion ratio, exit velocity, ambient pressure and exit 
pressure Since the STME uses a translatable skirt, the expansion ratio of the engine 
be set at'discrete values of35:1, 70:1, and 150:1. To insure that the engine* ope 
efficiently during all points in the launch, it was necessaiy to calculate the altitudes at 
which the skirt should be translated. This analysis was performed by calculating th 
coefficients for all three nozzles and finding the points of intersection. 

These calculations were based on the data presented in the Aeroject contractor's report for 
stamUrd atmosphere was used as a modal for pressure and dens.ty as 
functions of altitude. Most of the engine data was solved using simple, iterative 
techniques, isentropic relations, and graphical methods. 



2.1 Analysis 

Given an expansion ratio for a particular nozzle and an appropriate value for the 
coefficient of specific heats, the exit Mach number may be determined. 


Assuming y = 1.24 for a mixture of liquid oxygen and liquid hydrogen in a 6:1 mass ratio 
and the first (fixed) nozzle expansion ratio of 35:1, this isentropic relation holds: 

At _ A— j I +■ 1*1 e- 1- "1j z I t'H 

/X* **«-L J 

This function may be solved iteratively to obtain exit Mach numbers for the nozzles. This 
was solved using the included iterative, and the following were obtained as the exit Mach 
numbers: 


Nnasrie Expansion Ratio 
35:1 
70:1 
150:1 


Frit Mach Number 
4.369 
4.931 
5.581 


Next, the flow in a nozzle may be assumed to be very nearly isentropic, such that the 
chamber pressure may be related to the nozzle exit pressure through the isentropic relation 

Given a chamber pressure of 2385 psia under normal power levels and 3125 psia under 
emergency conditions, the following exit pressures were obtained: 


Pe (NFL) PeflSPLl 
Fore = 35:1 5.069 6.642 

For e = 70:1 2.058 2.696 

Fore =150:1 0.771 1.010 


For a given expansion ratio, a nozzle will have a maximum efficiency at a certain design 
altitude. At this altitude, the exit pressure is equal to the ambient pressure of the 
atmosphere. From the 1962 standard atmosphere, the design altitudes were determined to be 

hd (Normal P ower Level) (EmerggnCY POWCT Lgrel al 

For e = 35:1 26,700 ft 20,400 ft 

For e = 70:1 45,800 ft 40,200 ft 

For e = 150:1 66,300 ft 60,600 ft 


In order to obtain the proper altitudes at which the sections of the nozzle should be extended, 
the following procedure is followed. 


The thrust coefficient is defined as: 

Or = 


Pt a* 


Since thrust varies with altutde, this expression may be rewritten as. 

Ct -- 



w 


Pr 


2 


Ac 
A * 


8 



For example, consider the nozzle with an expansion ratio of 70:1 at an 


altitude of 50,000 feet. 


Cj ■■ * (j^f) 


11 355 


2-050 * 


- ( 10 ) 


CIS 


Similarly, thrust coefficients for other nozzle and atmospheric conditions may t>e 
calculated. As shown in Figure 2.1.1, the points where these graphs cross indicate the 
altitudes where the nozzles should be translated for maximum efficiency. 

For Normal Power Levels, the translation altitudes were 56 ’ 50 ° 

feet For Emergency Power Levels, the altitudes were 31,100 feet and 50,800 feet. 
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1000 


2000 


3000 


4000 


5000 


6000 


7000 


8000 


9000 


10000 


11000 


12000 


13000 


14000 


15000 


16000 


17000 


18000 


19000 


20000 


21000 


22000 


23000 


24000 


25000 


26000 


27000 


28000 


29000 


30000 


31000 


32000 


33000 


34000 


35000 


36000 


37000 


38000 


39000 


40000 


41000 


42000 


43000 


44000 


45000 


46000 


47000 


48000 


49000 


50000 


51000 


52000 


53000 


54000 


55000 


B 


P 


14.7000 


14 .1700 


13.6600 


13.1700 


12.6900 


12.2300 


11.7800 


11.3400 


10.9200 


10.5000 


10.1100 


9.7200 


9.3460 


8.9840 


8.6330 


8.2940 


7.9650 


7.6470 


7.3390 


7.0410 


6.7540 



UP h ,1*3 5 


1.6185 


1.62631 


1.6338 


1.64101 


1.6480 


1.6548 


1.6614 


1.6679 


1.67401 


1.6802 


1.685 91 


6.20701 


5.94701 


5.6960 


5.45401 


5.2200 


4.994 ol 


4.7770 


4.5670 


4.36401 


4.1690 


3,9810 


3,80001 


3.6260 


3.4580 


3.29701 


3.1420 


2.99401 


2.8540 


2.7200 


2.5920 


2.4710 


2.35501 


2.1390 


2 . 0390 l 


1.9430 


1.8520 


1.7650 


1.6820 


1.60301 


1.52801 


1.4560 


1.3880 


1.3230 


1.6971 


1.7024 


1.70761 


1.71261 


1.7174 


1.7220 


1.72661 


1.7309 


1.7352 


1.7392 


1.7432 


1.74701 


1.7507 


1.75421 


1.7577 


1.7610 


1.7642 


1.7672 


1.7702 


1.7731 


1.7758 


1.77851 


1.7811 


1.7835 


1.78591 


1.7882 


1.79031 


1.7924 


1.79431 


1.7962 


1.79801 


1.7997 


1.8013 


1.8029 


1.80431 


1.8058 


1.8071 


1.8084 


1.80961 


1.8107 


1.8118 


1.8129 


1.8139 


1.81491 


1 . 6696 


1.6756 


1.6813 


1.6868 


1.6921 


1.6973 


1.7023 


1.7073 


1.7120 


1.7167 


1.7210 


1.7254 


1.7296 


1.7336 


1.7376 


1.7414 


1.7451 


1.7486 


1.7521 


1.7554 


1.7586 


1.7617 


1.7647 


1.7677 


1.7705 


1.7732 


1.7758 


1.7783 


1.7808 


1.7831 


1.7854 


1.7876 


1.7897 


1.7917 


1.7937 


1.7955 


1.7973 


1.7991 


1.8007 


1.8023 


1.8038 


1.8052 


1.8066 


1.8079 


1.8091 


1.8103 


1.8114 


1.8125 


1.8135 


1.8145 


1.8154 


1.8163 


1.8172 


1.8180 


1.8187 


1.8194 


1.4649 


1.4798 


1.4942 


1.5083 


1.5218 


1.5350 



1.5840 


1.59551 


1.6064 


1.6171 


1.6274 


1.63731 


1.6470 


1.6563 


1.66531 


1.68251 


1.69071 


1.6986 


1.7062 


1.7136 


1.72071 


1.73421 


1.7405 


1.7467 


1.75271 


1.7584 


1.7639 


1.76921 


1.7793 


1.7840 


1.78851 


1.7929 


1.79701 


1.80091 


1.8047 


1.8082 


1.81161 


1.81491 


1.8180 


1.82091 


1.8237 


1.8264 


1.8289 


1.8314 


1.8337 


1.83591 


1.8380 


1.8400 


1.84191 


1.5515 


1.5633 


1.5748 


1.5857 


1 . 5965 L 


1.60681 


H 


£■150 


1.2164 


1.6267 


1.6361 


1.6455 


1.6543 


1.66301 


1.6795 


1.6874 


1.6950 


1.7023 


1.7094 


1.7163 


1 . 7295 ) 


1.73571 


1.7532 


1.7586 


1.7638 


1.76891 


1.7784 


1.7830 


1.7874 


1.7916 


1.7956 


1.7995 


1.8033 


1.8069 


1.8104 


1.8137 


1.8168 


1.8198 


1.8227 


1.8254 


1.8280 


1.8305 


1.8328 


1.8351 


1.8372 


1.8393 


1.8412 


1.8431 


1.8448 


1.8465 


1.8481 


1.8497 


1.8511 


1.0308 


1.0629 


1.0937 


1.1239 


1.1529 


1.1812 


1.2088 


1.2352 


1.2617 


1.2862 


1.3107 


1.3342 


1.3570 


1.3791 


1.4004 


1.4605 


1.4792 


1.4973 


1.5148 


1.5317 


1.54 80 


1.5638 


1.5790 


1.5937 


1.6079 


1.6216 


1.6348 


1.6476 


1.6598 


1.6717 


1.6830 


1.6940 


1.7046 


1.75101 


1.7590 


1.7666 


1 . 7809 ) 


1 .7875 


1 . 7938 


1.7998 


1.8056 
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Figure 2.1.1: Plot of Thust Coefficients vs. Altitude (NPL) 
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Figure 2.1.2: Thust Coefficients vs. Altitude (EPL) 
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3.0 Mission Analysis 

To prove the flight worthiness and technical feasibility of the ASTS as proposed, a launch 
profile was constructed and simulated. This simulation consisited of a dynamic model of 
the vehicle’s instantaneous position, velocity, acceleration, thrust, weight, orientation, 
and flight path with respect to the horizon. In addition to this vehicle-specific data, 
atmospheric conditions such as temperature, pressure, and density were modeled. 

The basic format of this simulation was a fourth-order Runge-Kutta integration scheme. A 
number of asumptions were necessary to develop this model of the launch trajectory. First, 
the model is constrained by a run time of approximate seven minutes and a total distance 
downrange not exceeding 800 nautical miles, so the surface of the earth was modeled as 
being flat. In addition, rotational effects of the earth were neglected. This assumption 
cannot degrade the validity of results, since it could only create conditions more favorable 
to a due east launch from Kennedy Space Center. Drag forces were neglected for two 
principal reasons: they can be considered small in comparison to thrust forces necessary 
for a ballistic trajectory, and they would not exceed the advantages gained by 
incorporating the earth's rotation. 

All initial conditions were set based on values calculated from the design process outlined 
in Section 1.0. Bum times, initial weights, total thrust available, and fuel available are the 
the most noteworthy examples. Thrust was computed as a funciton of altitude, based on 
translation altitudes for the extendable skirt on the STME nozzles, pressure forces, and 
known values of sea level and vacuum values. It was estimated that the maximum 
acceptable pitch rate was 0.5 degrees per second to retain stability during launch. This is 
the value used in the program for all pitch manuevers. 

Before finalizing the launch profile, it was necessary to establish the range of orbits into 
which the orbiter could be successfully launched. Based on current STS mission analyses 
and estimations of maximum altitude and velocity available at this altitude, our launch 
goal was projected to be an altitude of 100 nautical miles and a velocity equal to that of the 
local circular speed. Inclination was varied during flight both to achieve this orbit and to 
place the booster in an alignment practical for separation. 








program AE447Prog; 

($MC68020+} 

{SMC68881+} 

{SSETC Elems881:-true) 


Written by Del Johnson (in TML Pascal II vl.00 under MPW v2 .02)__ 

Written”for AE447-AE449 Shuttle II design group (Fall 19 88 - Spri ng 198in 


This program is designed to model the flight of design configurations for the 
Shuttle II Launch System. Design data was used in the model- *® S f. tune the 
execution of these programs was used in an iterative process to fine tune the 

model as presented in the report. 


Definition of variables u[i]: 


( Integral variables- } 


u [ 1] 
u( 2] 
u[ 3] 
u [ 4] 


u[ 
u ( 


5] 

6 ] 


- u[x] - Location along the x-axis 
» u(xd] — Time rate of change of location along the 
m u [y] — Location along the y-axis 

„ u [yd] ■ Time rate of change of location along the 
m u[th] — Angle with respect to the y-axis 
« u[thd] « Time rate of change of angle with respect 


x-axis 

y-axis 

to the y-axis 


( Other variables— } 
u( ] - u[alpha]- Angle of attack 

u[ ] - u[gamma]- Flight path angle 

u[ ] - u[wt] - Total vehicle weight 

u[ ] - u[wtd] - Time rate of change of total vehicle weight 

u( ] - u[dens] - Density of the atmosphere 

u[ ] - u[press] - Density of the atmosphere 
u[ ] - u[ thrust] - Thrust produced by all engines firing 


*) 

uses MemTypes,SANE; 


const 


{ Number of equations- } 

numlntEquations - 6; 

numTotEquations - 11; 


{ Time Considerations- ) 

tbegin - 0.00; [ seconds } 

tend - 420.00; { seconds } 

dt - 0.10; { seconds ] 

cint - 10? 

cols - 4; 

{ Subscripts of the variables of motion— } 


) 
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xd “ 1 * 

x ■ 2; 

yd “ 3; 

y - 4; 

wt “ 5; 

th - 6; 

vtot ” 7; 

gamma “ 8; 

thrust “ 9; 

press ” 10; 

dens ” 11; 

Separation “ 1 
Shutdown ” 2 


StartPitchl - 3; 

StopPitchl ■ 4; 

StartPitch2 - 5; 

StopPitch2 • 6; 

{ Constants related to the problem.- ) 

{ Mathematical constants.- } 

conv - 57.29577951; { degrees/radians } 


( 


Physical 

Re 

gO 

rho 

rhoO 

Psl 


constants.- } 

- 20938912.00000; ( feet } 

- 32.17400; { feet/second~2 } 

- 0.0023768; { pound-mass/feet*3 ) 

- 0.00237 68; { pound-mass/feet^ } 

- 14.70000; { pound-force/inches / '2 l 


{ Vehicle characteristics- ) 

thrO - 363535.00000; { lbf } 

wt d0 - -949.40000; { lbm/sec } 


pressAlt “ 25000 
densAlt “ 25000 

Ael - 3247 
Ae2 - 6503 
Ae3 - 14526 


00000; { feet } 
00000; { feet } 

22085; { inchest ) 
88219; { inches A 2 } 
72443; { inchest } 


designAlt1NPL 
designAlt2NPL 
designAlt3NPL 
designAlt1EPL 
designAlt2EPL 
designAlt3EPL 


26700.00000; { feet } 
45800.00000; { feet } 
66300.00000; { feet } 
20400.00000; { feet } 
40200.00000; { feet } 
60600.00000; { feet } 


transAltlNPL - 36900.0000; { feet } 
transAlt2NPL - 56500.0000; { feet } 
transAltlEPL - 31100.0000; { feet } 
transAlt2EPL - 50800.0000; ( feet } 


PelNPL - 5.069; { lbf/feet~2 } 
Pe2NPL - 2.058; { lbf/feet~2 } 
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Pe3NPL - 0.771; { lbf/feet A 2 } 
PelEPL - 6.642; { lbf/feet A 2 } 
Pe2EPL - 6.296; { lbf/feet A 2 } 
Pe3EPL - 1.010; { lbf/feet A 2 } 


thrMaxlNPL * 411259.15000; 
thrMax2NPL - 424577.04000; 
thrMax3NPL - 435000.00000; 

thrMaxlEPL * 545789.15000; 
thrMax2EPL - 556312.00000; 
thrMax3EPL - 580000.00000; 

{ Problem intial conditions... } 


xdO 

- 

0.00 

{ 

feet/second 

1 

xO 

• 

0.00 

{ 

feet } 


ydO 

- 

0.00 

{ 

feet/second 

} 

yO 

- 

0.00 

{ 

feet } 


thO 


90.00 

{ 

degrees } 


vtotO 

- 

0.00 

{ 

feet/second 

} 


type 

uArray 
name Array 
printrange 


array 
array 
1. .32; 


[1..numTotEquations] 
[1..numTotEquations] 


of extended; 
of str255; 


var 


outfile : text; 
tab : char; 

3 p : array [1..numTotEquations] of integer; 

DryWtBooster, DryWtOrbiter,wtO : extended; 
t : extended; 

time : array [1..6] of extended; 

Jc ; array [ 1 . . 4 ,1. .numlntEquations] of extended; 

u,uu : uArray; 

i, j,m,n : integer; 
printset : set of printrange; 


{ Variables of motion^ } 

wtd, thr, thd : extended; 

thrMax : extended; 


function g(alt : extended): extended; 
begin 

g g0*sqr (Re)/sqr (Re+alt) ; 

end; 

function density(alt : extended): extended; 
begin 

density :* rho0*exp(-alt/densAlt); 


end; 

function pressure(alt : extended): extended; 
begin 

pressure :« Psl*exp(-alt/pressAlt); 

end; 



17 








function pitchRate(t 
begin 


extended): extended; 


(t < Time[StartPitchl]) then pitchRate :« 
if (t < Time [StopPitchl] ) then pitchRate 

if (t < Time[StartPitch2]) then pitchRate 

if (t < Time[StopPitch2] ) then pitchRate 

if (t < tend) then pitchRate 0.00; 


end; 

function thrustAvail(alt 


extended): extended; 


0.00 else 

-0.50 else 
:« 0.00 else 

-0.50 else 


var 

Ae,thrl : extended; 

begin 

if (alt<transAltlNPL) then begin 
Ae Ael; 
thrl :» thrMaxlNPL; 


end 

else 

if (alt<transAlt2NPL) then begin 
Ae :■ Ae2; 
thrl :» thrMax2NPL; 

end 

else begin 

Ae :• Ae3; 

thrl :« thrMax3NPL; 

end; 

thrustAvail :* thrl - pressure(alt)*Ae; 

end; 

(———-——————----—-~i 

procedure initVariables; 
begin 

t :■ tbegin; 
u[yd] ydO; 
u[y] yO; 

u[xd] :■ xdO; 
u[x] xO? 

u[wt] :* wtO; 
u[th] :■ thO; 
u{gamma] 90; 

u[vtot] sqrt(sqr(u[yd])+sqr(u(xd])); 
u[press] pressure(u[y]); 

u(dens] :« density(u[y]); 

end; 

procedure print(style ; integer); 
var 

i : integer; 
begin 

case style of 
1 : begin 

write(outfile,'t - ',t:6:2); 

for i :* 1 to numTotEquations do begin 

if (i in printset) then write(outfile,' u[',i:2,'] 

,u[i]:12:2) ; 

if (i in printset) then write( ' u[',i:2,'] 

-',u[i]:12:2); 

end; 
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') ; 
'); 


writeln{outfile); 
writeln; 
end; 

2 ; begin 

i 0; 

writeln(outfile,'Time - ',t:7:2); 
writeln( 'Time m ’,t:7:2); 

for j :■ 1 to numTotEquations do 
if (j in printset) then begin 
i i+1; 

write (outfile,'u[\ j:2,'] ,u[j]:12:2,' 

write( 'u[',j:2,'] “',u[j]:12:2,' 

if (i mod cols - 0) then begin 
i :« 0; 

writeln(outfile); 
writeln; 
end; 
end; 

if (i<>0) then begin 
writeln(outfile); 
writeln; 
end; 
end; 

3 : begin 

write(outfile,t:7:2,tab); 

for j :■ 1 to numTotEquations do 

if (j in printSet) then write (outf ile,u[j]:12:sp[j], tab); 
writeln(outfile); 
end; 
otherwise; 
end; (case) 

end; mm mm , _ 

{ Total function defining variables in terms of one another.™ ) 
function f(j: integer; t: extended; u: uArray): extended; 
begin 


case j of 
xd 
x 

yd 

y 

wt 
th 


f 
f 
f 
f 
f 
f 

otherwise 


- thr*cos(u[th]/conv)*g(u[y])/u[wt]; 

- u[xd]; 

- thr*sin(u[th]/conv)*g(u[y])/u[wt] - g(u[y]); 

- u[yd]; 

* wtd; 

thd; 


end; 


end; 


begin 


Textbook(nil); 
tab ;" chr(09); 

for j :■ 1 to numTotEquations do 
sp[j] :- 2; 
sp[press] :■ 8; 
sptdens] :■ 8; 


open (outfile,'l/AE448prog.dat ); 
rewrite(outfile); 
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m numlntEquations; t Number of integration equations } 

n trunc((tend-tbegin)/dt); { Number of steps } 

{ Set names of all variables.- } 

printset [xd,x,yd,y, wt,vtot,th,gamma,thrust,press,dens]; 

{ Set design variables— } 

wtO 2769714.29; { lbf } 

DryWtBooster :~ 170000.00; { lbf } 

DryWtOrbiter 260000.00; { lbf } 

Time[Separation] :■ 148.00; { seconds } 

Time[Shutdown] :« 394.00; { seconds } 

Time[StartPitchl] 20.00; { seconds } 

Time[StopPitchl] 140.00; ( seconds } 

Time[StartPitch2] 160.00; { seconds } 

Time[StopPitch2] :■ Time[StartPitch2] + 38.00; { seconds } 

[ Set initial conditions of all variables- } 
initVariables; 

[ Echo initial conditions of all variables in the printing set... > 
print(3); 

{ Beginning of time-step loop— ) 
for i 1 to n do begin 

{ Computation of other, non time-dependent variables- } 

if (t ■ Time[Separation]) then u[wt] :* u[wt] - DryWtBooster 
thrMax thrustAvail(u[y]); 

if (t < Time[Separation]) then thr ;■ 10*thrMax else 
if (t < Time[Shutdown] ) then thr 4*thrMax else 

if (t < tend) then thr 0*thrMax; 

u[press] :« pressure(u[y]); 
ujdens] density(u[y]); 

u[thrust] ;■ thr; 

if (t < Time(Separation]) then wtd :■ 10*wtd0 else 

if (t < Time[Shutdown] ) then wtd 4*wtd0 else 

if (t < tend) then wtd 0; 

thd pitchRate(t); 

if (u[th] > 360.00) then u[th] u[th] - 360.00; 
if (u[thj < -360.00) then u[th] :« u(th] + 360.00; 

( Runge-Kutta 4th-Order integration loop of all variables- ) 

for j :■ 1 to m do 

k(1,j] dt*f(j,t,u); 
for j :■ 1 to n do 

uu(j] :“ u[j] + k[l,j]/2; 
for j 1 to m do 

k[2,j] dt * f(j,t+dt/2, uu); 
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for j 1 to m do 

uu[j] u[j] + k[2,j]/2; 

for j :» 1 to in do 

k[3,j] dt*f(j,t+dt/2,uu); 

for j 1 to in do 

uu[j] : m u(j] + k13 , j]; 

for j 1 to m do 

k[4,j] dt*f(j,t+dt,uu); 

for i 1 to m do . 

u [j] u[j] + {k[l,j]+2*k[2,j]+2*k[3,Dl+k[4 # Dl)/6; 

t :■ tbegin + i*dt; 

u[vtot] :■ sqrt(sqr (u[xd])+sqr(u[yd])); 

if (u[xd] > 0) then u[gamma] :■ arctan(u[yd]/u[xd]) *conv 
u[gamma] :■ 90.00; 

{ Printing parameters™ ] 

if (i mod cint - 0) and (printset <> []) then print(3); 


end; 

close(outfile); 


end. 


else 
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l/AE 448 data 



105.00 


106.00 


107.00 


108.00 


109.00 


110.00 


111.00 


112.00 


113.00 


114.00 


115.00 


116.00 


117.00 


118.00 


119.00 


120.00 


121.00 


122.00 


123.00 


124.001 


125.00 


126.00 


127.00 


128.00 


129.00 


130.001 


131.00 


132.00 


133.00 


134.00 


135.00 


136.00 


137.00 


138.00 


139.00 


140.00 


141.00 


142.00 


143.00 


144.00 


145.00 


146.00 


147.00 


148.00 


149.00 


150.00 


151.001 


152.00 


153.00 


154.001 


155.00 


156.00 


157.00 


B 


1940.72 


1993.34 


2046.78 


2101.05 


2156.14 


2212.06 


2268.81 


2326.39 


2384.81 


2444.07 


2504.16 


2565.10 


2626.88 


2689.52 


2753.00 


2817.33 


2882.53 


2948.58 


3015.50 


3083.28 


3151.94 


3221.47 


3291.88 


3363.17 


3435.34 


3508.41 


3582.37 


3657.23 


3733.00 


3809.67 


3887.26 


3965.76 


4045.19 


4125.54 


4206.83 


4289.05 


4372.01 


4455.50 


4539.54 


4624.12 


4709.25 


4794.94 


4881.21 


4968.05 


5007.92 


5047.90 


5088.01 


5128.23 


5168.57 


5209.03 


5249.62 


5290.32 


5331.15 


C 


51134.15 


53101.11 


55121.10 


57194.95 


59323.47 


61507.50 


63747.86 


66045.40 


68400.93 


70815.30 


73289.34 


75823.90 


78419.83 


81077.96 


83799.14 


86584.24 


89434.10 


92349.58 


95331.55 


98380.86 


101498.40 


104685.03 


107941.63 


111269.08 


114668.26 


118140.06 


121685.38 


125305.10 


129000.14 


132771.40 


136619.79 


140546.22 


144551.61 


148636.90 


152803.00 


157050.86 


161381.35 


165795.06 


170292.53 


174874.31 


179540.95 


184293.00 


189131.03 


194055.61 


199043.58 


204071.48 


209139.43 


214247.54 


219395.93 


224584.72 


229814.04 


235084.00 


240394.72 


2049.44 


2074.47 


2099.40 


2124.20 


2148.89 


2173.431 


2197.83 


2222.07 


2246.15 


2270.06 


2293.78 


2317.30 


2340.63 


2363.74 


2386.64 


2409.31 


2431.74 


2453.93 


2475.86 


2497.53 


2518.93 


2540.05 


2560.88 


2581.42 


2601.65 


2621.57 


2641.16 


2660.43 


2679.35 


2697.93 


2716.16 


2734.02 


2751.50 


2768.61 


2785.32 


2801.63 


2817.90 


2834.48 


2851.38 


2868.61 


2886.17 


2904.05 


2922.28 


2940.84 


2932.30 


2923.83 


2915.44 


2907.13 


2898.90 


2890.74 


2882.66 


2874.66 


2866.74 


91838.63 


93900.60 


95987.54 


98099.35 


100235.90 


102397.07 


104582.72 


106792.68 


109026.81 


111284.93 


113566.86 


115872.41 


118201.401 


120553.60 


122928.81 


125326.811 


127747.35 


130190.21 


132655.121 


135141.84 


137650.09 


140179.60 


142730.09 


145301.26 


147892.82 


150504.451 


153135.84 


155786.66 


161145.26 


163852.33 


166577.451 


169320.24 


172080.33 


174857.331 


177650.84 


180460.57 


183286.741 


186129.64 


188989.61 


191866.97 


194762.05 


197675.19 


200606.721 


203543.29 


206471.35 


209390.98 


212302.26 


215205.27 


218100.08 


220986.77 


223865.42 


226736.11 


1772844.29 


1763350.29 


1753856.29 


1744362.29 


1734868.29 


1725374.29 


1715880.29 


1706386.29 


1696892.29 


1687398.29 


1677904.29 


1668410.29 


1658916.29 


1649422.29 


1639928.29 


1630434.29 


1620940.29 


1611446.29 


1601952.29 


1592458.29 


1582964.29 


1573470.29 


1563976.29 


1554482.29 


1544988.29 


1535494.29 


1526000.29 


1516506.29 


1507012.29 


1497518.29 


1488024.29 


1478530.29 


1469036.29 


1459542.29 


1450048.29 


1440554.29 


1431060.29 


1421566.29 


1412072.29 


1402578.29 


1393084.29 


1383590.29 


1374096.29 


1364602.29 


1190804.69 


1187007.09 


1183209.49 


1179411.89 


1175614.29 


1171816.69 


1168019.09 


1164221.49 


1160423.89 
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l / AE 448 data 



158.00 


159.00 


160.00 


161.00 


162.00 


163.00 


164.00 


165.00 


166.00 


167.00 


168.00 


169.00 


170.00 


171.00 


172.00 


173.00 


174.00 


175.00 


176.00 


177.00 


178.00 


179.00 


180.00 


181.00 


182.00 


183.00 


184.00 


185.00 


186.00 


187.00 


188.00 


189.00 


190.00 


191.00 


192.00 


193.00 


194.00 


195.00 


196.00 


197.00 


198.00 


199.00 


200.00 


201.00 


202.00 


203.00 


204.00 


205.00 


206.00 


207.00 


208.00 


209.00 


210.00 


B 


5372.10 


5413.17 


5495.79 


5537.55 


5579.64 


5622.06 


5664.81 


5707.89 


5751.29 


5795.02 


5839.08 


5883.45 


5928.15 


5973.17 


6018.51 


6064.17 


6110.14 


6156.43 


6203.04 


6249.95 


6297.18 


6344.72 


6392.57 


6440.72 


6489.18 


6537.94 


6587.01 


6636.37 


6686.04 


6736.00 


6786.26 


6836.81 


6887.66 


6990.23 


7041.94 


7093.94 


7146.23 


7198.79 


7251.64 


7304.72 


7357.99 


7411.46 


7465.11 


7518.97 


7573.01 


7627.26 


7681.71 


7736.36 


7791.21 


7846.27 


7901.53 


245746.33 


251138.95 


256572.71 


262047.76 


267564.41 


273122.98 


278723.80 


284367.21 


290053.53 


295783.09 


301556.22 


307373.25 


313234.48 


319140.26 


325090.89 


331086.71 


337128.02 


343215.15 


349348.41 


355528.12 


361754.59 


368028.13 


374349.06 


380717.67 


387134.29 


393599.21 


400112.75 


406675.19 


413286.86 


419948.04 


426659.03 


433420.14 


440231.65 


447093.87 


454007.07 


460971.56 


467987.62 


475055.53 


482175.59 


489348.08 


496573.27 


503851.44 


511182.78 


518567.49 


526005.76 


533497.78 


541043.75 


548643.88 


556298.34 


564007.36 


571771.13 


579589.85 


587463.73 


2858.891 


2851.13 


2843.45 


2835.671 


2827.60 


2819.25 


2810.61 


2801.681 


2792.45 


2782.92 


2773.081 


2762.94 


2752.48 


2741.71 


2730.63 


2719.22 


2707.48 


2695.42 


2683.02 


2670.28 


2657.21 


2643.79 


2630.02 


2615.90 


2601.43 


2586.60 


2571.41 


2555.85 


2539.92 


2523.61 


2506.93 


2489.87 


2472.42 


2454.59 


2436.36 


2417.74 


2398.72 


2379.29 


2359.46 


2339.21 


2297.70 


2276.90 


2256.14 


2235.421 


2214.75 


2194.13 


2173.551 


2153.01 


2132.52 


2112.07 


2091.68 


2071.32 



229598.92 


232453.93 


235301.21 


238140.79 


243795.90 


246610.86 


249417.03 


252214.12 


255001.83 


257779.85 


260547.89 


266052.75 


268788.95 


271513.89 


279617.99 


282294.67 


284958.44 


287608.97 


290245.91 


292868.90 


295477.60 


298071.64 


300650.68 


315788.96 


318252.50 


320698.01 


325 S 33.36 


327922.40 


330291.80 


332641.17 


334970.09 


344077.80 


346302.89 


348507.33 


350691.16 


359221.36 


361302.85 
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l / AE 440 data 



211.00 


212.00 


213.00 


214.00 


215.00 


216.00 


217.00 


218.00 


219.00 


220.00 


221.00 


222.00 


223.00 


224.00 


225.00 


226.00 


227.00 


228.00 


229.00 


230.00 


231.00 


232.00 


233.00 


234.00 


235.00 


236.00 


237.00 


238.00 


239.00 


240.00 


241.00 


242.00 


243.00 


244.00 


245.00 


246.00 


247.00 


248.00 


249.00 


250.00 


251.00 


252.00 


253.00 


254.00 


255.00 


256.00 


257.00 


258.00 


259.00 


260.00 


261.00 


262.00 


263.00 


7957.00 


8012.68 


8068.58 


8124.69 


8181.01 


8237.55 


8294.31 


8351.29 


8408.50 


8465.93 


8523.59 


8581.47 


8639.59 


8697.94 


8756.53 


8815.35 


8874.41 


8933.72 


8993.27 


9053.06 


9113.10 


9173.40 


9233.94 


9294.74 


9355.80 


9417.12 


9478.70 


9540.54 


9602.66 


9665.04 


9727.69 


9790.62 


9853.83 


9917.31 


9981.08 


10045.14 


10109.48 


10174.11 


10239.04 


10304.27 


10369.79 


10435.62 


10501.75 


10568.19 


10634.94 


10702.01 


10769.39 


10837.10 


10905.13 


10973.49 


11042.18 


11111.20 


11180.56 


C 


595392.98 


603377.80 


611418.41 


619515.03 


627667.86 


635877.12 


644143.03 


652465.81 


660845.69 


669282.88 


677777.62 


686330.13 


694940.64 


703609.39 


712336.60 


721122.52 


729967.38 


738871.42 


747834.89 


756858.03 


765941.09 


775004.32 


784287.97 


793552.29 


802877.54 


812263.98 


821711.86 


831221.46 


840793.04 


850426.86 


860123.21 


869882.34 


879704.54 


889590.09 


899539.26 


909552.35 


919629.64 


929771.41 


939977.96 


950249.59 


960586.59 


970989.27 


981457.93 


991992.87 


1002594.41 


1013262.86 


1023998.53 


1034801.75 


1045672.84 


1056612.12 


1067619.92 


1078696.53 


2051.02 


2030.76 


2010.541 


1990.38 


1970.26 


1950.19 


1930.17 


1910.20 


1890.28 


1870.40 


1850.58 


1830.81 


1811.08 


1791.41 


1771.79 


1752.22 


1732.70 


1713.23 


1693.81 


1674.45 


1655.14 


1635.89 


1616.68 


1597.54 


1578.44 


1559.41 


1540.42 


1521.50 


1502.63 


1483.81 


1465.06 


1446.36 


1427.72 


1409.14 


1390.61 


1372.15 


1353.75 


1335.40 


1317.12 


1298.90 


1280.74 


1262.65 


1244.62 


1226.65 


1208.74 


1190.90 


1173.13 


1155.42 


1137.77 


1120.20 


1102.69 


1085.25 


1067.88 


363364.02 


365404.90 


367425.55 


369426.01 


371406.32 


373366.55 


375306.73 


377226.91 


379127.15 


381007.49 


382867.97 


384708.66 


386529.60 


388330.85 


390112.44 


391874.44 


393616.89 


395339.85 


397043.36 


398727.49 


400392.28 


402037.79 


403664.07 


405271.18 


406859.16 


408428.08 


409977.99 


411508.95 


413021.01 


414514.22 


415988.65 


417444.36 


418881.39 


420299.81 


421699.68 


423081.06 


424444.00 


425788.57 


427114.83 


428422.84 


429712.66 


430984.35 


432237.98 


433473.60 


434691.29 


435891.11 


437073.12 


438237.38 




440512.95 


441624.391 


442718.351 


951555.89 


947758.29 


943960.69 


936365.49 


932567.89 


928770.29 


924972.69 


921175.09 


913579.89 


909782.29 


905984.69 


902187.09 


898389.49 


894591.89 


890794.29 


886996.69 


883199.09 


879401.49 


875603.89 


871806.29 


868008.69 


864211.09 


860413.49 


056615.89 


852818.29 


849020.69 


845223.09 


841425.49 


837627.89 


833830.29 


830032.69 


826235.09 


822437.49 


818639.89 


814842.29 


811044.69 


807247.09 


803449.49 


799651.89 


795854.29 


792056.69 


788259.09 


784461.49 


780663.89 


776866.29 


773068.69 


769271.09 


765473.49 


761675.89 
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264.00 


265.00 


266.00 


267.00 


268.00 


269.00 


270.00 


271.00 


272.00 


273.00 


274.00 


275.00 


276.00 


277.00 


278.00 


279.00 


280.00 


281.00 


282.00 


283.00 


284.00 


285.00 


286.00 


287.00 


288.00 


289.00 


290.00 


291.00 


292.00 


293.00 


294.00 


295.00 


296.00 


297.00 


298.00 


299.00 


300.00 


301.00 


302.00 


11320 . 


11390 . 


11461 . 


11532 . 


11604 . 


11675 . 


11748 . 


11820 . 


11893 . 


11966 . 


12040 . 


12114 . 


12189 . 


12264 . 


12339 . 


12415 


12491 


12567 


12644 


12721 


12799 


12877 


12956 


13035 


13115 


13195 


304.001 


305.00 


306.00 


307.00 


308.00 


310.00 


311.00 


13356 


13438 


13520 


13602 


13685 


313.00 


314.001 


315.00 


316.001 


13852 . 


13936 . 


14021 . 


14107 . 


14192 . 


14279 . 


14366 . 


14453 . 


14542 , 


14630 


14719 


14809 


14900 


14991 


15082 


15174 


15267 


15361 


15455 


1101057.82 


3111112343.08 


. 71111236 98.56 


.46 


.571 1146 621.60 


.03 


. B 6 i 1169829.78 


. 061 1181541.71 


i - 62 l 1193326.02 


l . 57 l 120 5183.08 


1.89 


I 601 1229116.99 


1.691 124 1194.60 


>.18 


i . 061 1265573.09 


>- 351 1277874.76 


s 04 l 1290251.92 


l . 14 i 1302704.98 


7.661 1315234.34 


4-591 1327840.43 


1.951 1340523.66 


9-741 1353284.47 


7.961 1366123.29 


6.621 137904 0.54 


5.73 


5.291 1405112.15 


5.301 1418267.40 


1431502.90 


6 . 7011444819.09 


ft .111 1458216.46 


f > OOl 1471695.48 


i 2.371 1485256.62 


15.221 1498900.37 


1512627.23 


\2 - 431 1526437.69 


16.791 1540332.25 


>1 .661 1554311.43 


) 7.051 15683 75.74 


1050.581 444854 .13 


1033.34 


446920.84 


999 091 447928.48 


982.081 448919.06 


965.131 449892.65 


450849.34 


931.461 451789.20 


qi a 7i 452712.30 




65 _ 

8 4 l 1743985.131 


,63 1759206.31 
05 1774520 . 60 l 


111 1789928.63 


865.03 


848 . 62 ] 


832.28 


816.02 


799.85 


783.75 


767.74 


751.81 


735.96 


456238.60 


457079.04 


457903.18 


458711.11 


460278.64 


461038.41 


461782.28 



462510.36 


463222.71 

688.93 

463919.43 

673.43 



642.69 


627.45 


612.31 


597.26 


582.30 


567.44 


552.67 


537.99 


523.42 


508.951 


494.57 


480.30 


466.13 


452.06 


1596761.86 

438.09 

473459.76 

1611084.74 

424.24 

473890.92 

1625494.88 

410.49 

474308.27 

1639992.83 

396.85 

474711.93 

1654579.16 

383.32 

475102.01 

1669254.43 

369.90l 

475478.61 

1684019.21 

356.59 

475841.85 

1698874.07 

343.40 

476191.83 

1713819.62 

330.33 

476528.69 


317.371 


304.531 


291.81 


279 . 22L 


266.751 



465266.31 


465916.66 


466551.72 


467171.59 


467776.37 


468366.14 


468941.00 


469501.04 


470046.36 


470577.06 


471093.24 


471594.99 


472082.41 


472 SS 5.61 


473014.70 


477163.47 


477461.64 


477747.14 


478020.11 
























































































































































































































































































































00 o\ 

04 04 
CO CO 


l / AE 448 data 



317.00 


318.00 


319.001 


320.00 


321.00 


322.00 


323.00 


324.00 


325.00 


326.00 


327.00 


328.00 


329.00 


330.00 


331.00 


332.00 


333.00 


334.00 


335 . 

,00 

336 . 

,00 

337 . 

.00 

338 . 

.00 

339 

.00 

340 

o 

o 

341 

.00 

342 

.00 


343.00 


344.00 


345.00 


346.00 


347.00 


348.00 


349.00 


350.00 


351.00 


352.00 


353.00 


354.00 


355.00 


356.00 


357.00 


358.00 


359.00 


360.00 


361.00 


362.00 


363.00 


364.00 


365.00 


366.00 


367.00 


368.00 


369.00 


B 


15549 


15645 


15741 


15837 


15935 


16033 


16131 


16231 


16331 


16432 


16534 


16636 


16739 


16843 


16948 


17053 


17160 


17267 


17484 . 


17593 . 


17704 . 


17815 . 


17928 . 


18041 . 


18155 . 


18270 . 


18386 . 


18504 . 


18622 . 


18741 . 


18861 . 


18982 . 


19105 . 


19228 . 


19352 , 


19478 , 


19605 


19733 


19862 


19992 


20124 


20257 


20391 


20527 


20663 


20802 


20941 


21082 


21225 


21369 


21514 


21661 


801 1805431.03 


151 1821028.45 


15 1836721.54 


fl 3 l 1852510.98 


18 1868397.42 


22 1884381.57 


96 l 1900464.10 


4 ll 1916645.72 


57 1932927.15 


47 l 1949309.11 


.10 1965792.33 


49 l 1982377.57 


.641 1999065.57 


. 561 2015857.11 


. 2712032752.96 


. 781 2049753.92 


.101 2066860.79 


. 2412084074.38 


2101395.54 


.041 2118825.091 


.731 2136363.90 


- 29 l 2154012.84 


.75 2171772.781 


- 1112189644.63 


.39 2207629.31 


.61 2225727.731 


.78 2243940.84 


. 9 ll 2262269.60 


.031 2280714.991 


.15 2299278.00 


- 29 l 2317959.631 


.351 2551785.16 


.53 2572109.49 


.031 2592568.66 


.87 2613163.99 


.07 2633896.84 


.68 2654768.60 


.71 2675780.68 


.20 2696934.51 


.171 2718231.56 


.65 2739673.35 


.691 2761261.39 


254.40 


242.18 


230.09 


218.13 


206.30 


194.61 


183.05 


171.64 


160.36 


149.22 


138.23 


127.38 


116.69 


106.14 


95.75 


85.51 


75.43 


65.51 


55.75 


46.16 


36.74 


27.49 


18.41 


9.50 


0.78 


- 7.77 


- 16.12 


- 24.29 


- 32.27 


- 40.06 


- 47.65 


2336760.93 

- 55.03 

2355682.92 

- 62.21 

2374726.68 

- 69.18 

2393893.28 

- 75.94 

2413183.83 

- 82.49 

2432599.45 

- 88.81 

2452141.291 

- 94.91 

2471810.50 

- 100.79 

2491608.27 

- 106.43 

2511535.82 

- 111.83 

2531594.36 

- 116.99 


478280.681 


478528.961 


478765.08 


478989.18 


479201.39 


479401.83 


479590.65 


479767.98 


480088.74 


480232.46 


480365.25 


480487.28 


480598.68 


480699.61 


480790.23 


480870.68 


480941.14 


481001.76 


481052.71 


481094.14 


481126.24 




481168.24 


481164.73 


481132.54 


481104.24 


481068.06 


481024.19 


480972.83 


480914.20 


480848.48 




480696.66 


480610.99 


480421.24 


480317.61 


552807.89 


549010.29 


545212.69 


- 121.91 


- 126.58 


- 130.991 


- 135.14 


- 139.02 


- 142.64 


- 145.97 


- 149.031 


- 151.79 


- 154.27 


- 156.44 


480094.03 


479974.56 


479850.29 


479721.49 


479588.40 


479451.30 


479310.45 


479166.12 


479018.59 


478868.16 


478715.10 


478559.73 
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370.00 


371.001 


372.00 


373.00 


374.00 


375.00 


376.00 


377.001 


378.00 


379.00 


380.00 


381.00 


382.00 


383.00 


384.00 


385.00 


386.001 


387.00 


388.00 


389.00 


390.00 


391.00 


392.001 


393.00 


394.00 


395.00 


396.00 


397.00 


398.00 


399.00 


400.00 


401.00 


402.00 


403.00 


404.001 


405.00 


406.00 


407.001 


408.00 


409.00 


410.001 


411.00 


412.00 


413.00 


414.00 


415.00 


416.001 


417.00 


418.00 


419.00 


420.00 


B 


21810. 


21960. 


22112. 


22266. 


22421. 


22578. 


22737. 


22897. 


23060. 


23225. 


23391. 


23560. 


23730. 


23903. 


24078. 


24255. 


24435. 


24617. 


24801. 


24988. 


25177. 


25369. 


25564. 


25761. 


25962. 


25962. 


25962. 


25962. 


25962. 


25962. 


25962. 


25962 


25962 


25962 


25962 


25962 


25962 
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25962 
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311 2782997.25 


54l 2804882.54 


43 2826918.89 


OOl 2849107.96 


3ll 2871451.47 


37 2893951.17 


25 2916608.82 


97 2939426.28 


59 2962405.40 


14l 2985548.10 


67 3008856.33 


241 3032332.12 
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-158.30 

478402.33 

-159.85 

478243.23 

-161.08 

478082.74 

-161.98 

477921.18 

-162.55 

477758.89 

-162.77 

477596.20 

-162.65 

477433.46 

-162.16 

477271.02 

-161.31 

477109.26 

-160.08 

476948.54 

-158.46 

476789.23 

-156.46 

476631.74 

-154.04 

476476.45 

-151.22 

476323.79 

-147.96 

476174.16 

-144.28 


, -140.14 

475885.76 

1 -135.55 

475747.87 

1 -130.48 

475614.82 

1 -124.94 

475487.07 

1 -118.89 

475365.11 

i -112.33 

475249.45 

! -105.25 

475140.62 

! -97.63 

475039.13 

1 -89.46 

474945.54 

l -120.22 

474840.70 

) -150.98 

474705.10 


.161 3431662.26 


.161 3457624.42 


.16 3483586.58 


.161 3509548.74 


.16 3535510.91 


.16 3561473.07 


.161 3587435.23 


.16 3613397.39 


.16 3639359.55 


.161 3665321.71 


.16 3691283.87 


.1613717246.03 


.16 3743208.19 


.161 3769170.36 


.16 3795132.52 


.161 3821094.68 


:.16 3847056.84 


1.161 3873019.00 


:.16 3898981.16 


! ,16 3924943.32 


1.16 3950905.48 


i.16 3976867.64 


!.16 4002829.81 


!.16l 4028791.97 


-181.74 


-212.51 


-243.27 


-274.04 


-335.57 


-366.34 


-397.11 


-427.88 


-458.65 


-489.42 


-520.20 


-550.97 


-581.75 


-612.53 


-643.31 


-674.09 


-704.88 


-735.66 


-766.45 


-797.24 


-828.04 


-858.83 


-889.63 


474538.74 


474341.61 


474113.72 


473855.07 


473565.64 


473245.46 


472894.50 


472512.78 


472100.28 


471657.02 


471182.98 


470678.17 


470142.58 


469576.22 


468979.08 


468351.16 


467692.46 


467002.98 


466282.71 


465531.65 


464749.80 


463937.16 


463093.73 


462219.50 


351535.09 


347737.49 


343939.89 


340142.29 


336344.69 


332547.09 


328749.49 


324951.89 


321154.29 


317356.69 


313559.09 


309761.49 


305963.89 


302166.29 


298368.69 


290773.49 


286975.89 


283178.29 


279380.69 


275583.09 


271785.49 


267987.89 
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260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 
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90.001 


90.00 


90.001 


90.00 


90.00 


90.001 


90.00 


90.001 


90.00 


90.001 


90.00 


90.00 


90.001 


90.00 


90.00 


90.001 


90.00 


90.00 


90.00 


90.001 


90.00 


89.50 


89.00 


88.50 


88.00 


87.50 


87.00 


86.501 


86.00 


85.50 


85.00 


84.50 


84.001 


83.50 


83.00 


82.501 


82.00 


81.50 


81.001 


80.50 


80.00 


79.501 


79.00 


78.50 


78.00 


77.50 


77.00 


76.501 


76.00 


75.50 


75.00 


74.50 


H 


Vtotall 


0.00 


10.13 


20.40 


30.831 


41.41 


52.15 


63.051 


74.11 


85.35 


96.751 


108.34 


120.10 


132.041 


144.18 


156.50 


169.02 


181.74 


194.66 


207.791 


221.13 


234.69 


248.461 


262.45 


276.67 


291.121 


305.80 


320.72 


335.891 


351.31 


366.99 


382.93 


399.15 


415.66 


432.451 


449.54 


466.94 


484.651 


502.69 


521.05 


539.751 


558.80 


578.191 


597.951 


618.07 


638.56 


659.431 


680.69 


702.34 


724.391 


746.84 


769.70 


792.97 
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90.00 


90.00 


90.00 


90.00 


90.00 


90.00 


90.00 


90.00 


90.00 


90.00 


90.00 


90.00 


90.00 


90.00 


90.00 


90.00 


90.00 


90.00 


90.00 


90.00 


90.00 


89.95 


89.82 


89.62 


89.36 


89.05 


88.69 


88.29 


87.86 


87.39 


86.91 


86.40 


85.87 


85.32 


84.76 


84.19 


83.60 


83.01 


82.41 


81.80 


81.19 


80.58 


79.96 


79.34 


78.72 


78.09 


77.47 


76.85 


76.22 


75.60 


74.98 


74.37 


3698357.601 


3703490.101 


12.7565671 


12.5985091 


3770407.26 


3777262.78 


3784212.51 


3805537.43 


3812773.65 


3820058.15 


3827381.39 


3834733.78 


10.537757 


10.326637 


10.112616 


9.895967 



8.556786 


8.329764 


8.102438 


0.002063 


0.002037 


0.002011 




0.001529 


0.001493 


0.001457 


0.001420 


0.001384 
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816.67 


73.50 

840.79 


73.00 

865.34 


72.50 

890.331 

71.92 

72.00 

915.75 

71.32 

71.501 

941.62 

70.72 
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967.94 


70.50 

994.71 


70.00 

1021.94 



1049.63 

68.36 


69.00 


68.001 


67.501 


67.00 


66.501 


66.00 


65.50 


65.001 


64.50 


64.00 


63.50 


63.001 


62.50 


62.001 
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61.00 


60.50 


60.001 


59.50 


59.00 


58.50 


58.00 


57.50 


57. 00L 


56.501 


56.00 


55.501 


55.00 


54.501 


54.001 


53.50 


53.00 


52.501 


52.00 


51.501 


51.00 


50.50 


50.001 


49.50 


49.00 


48.50 


48.001 


1077.78 


1106.40 


1135.49 


1165.051 


1195.09 


1225.61 


1256.61 


1288.09 


1320.06 


1352.53 


1385.80 


1419.72 


1454.21 


1489.28 


1524.92 


1561.151 


1597.95 


1635.32 


1673.28 


1711.81 


1750.91 


1830.86 


1871.70 


1912.951 


1954.76 


1997.27 


2040.47 


2084.35 


2128.92 


2174.161 


2220.07 


2266.66 


2313.91 


2361.82 


2410.40 


2459.631 


2509.52 


2560.06 


2611.251 


2663.10 


2715.59 


2768.731 


67.77 


67.20 


66.62 


66.06 


65.49 
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62.19 


61.65 


6 


60.60 


60.07 


59.56 


59.04 


58.53 


58.03 


57.52 


56.53 


56.04 


55.56 


55.08 


54.60 


54.12 


53.65 


53.18 


52.72 


52.26 


51.80 


51.35 


50.90 


50.45 


50.01 


49.56 


49.13 


48.69 


48.26 


47.83 


47.40 


46.98 
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2822.52 


2876.95 


2987.75 


3044.12 


3101.13 


3158.79 


3217.10 


3335.66 


3456.82 


3580.61 


3643.49 


3707.04 


3771.25 


3836.13 


38.00 


34.50 


3 


33.50 


32.50 


3967.91 


4034.81 


4102.40 


4170.68 


4239.65 


4309.31 


4379.67 


4450.74 


4522.52 


4595.02 


4668.24 


4742.18 


31.00 


30.50 


3 


5045.33 


5123.00 


FTF1I 

30.00 

5201.44 

EEQI 

30.00 

5280.70 


3 


30.00 


30.001 


3 


30.00 


30.00 


30.00 


30. OOl 


30.00 


30.001 


3 


30 . OOl 


30.00 


30. OOl 


30.00 


5360.76 


5441.63 


5523.31 


5605.80 


5689.10 


5773.22 


5803.24 


5833.54 


5864.10 


5894.93 


5926.02 


5957.38 


5989.01 


6020.89 


6053.041 


40.54 


37.14 


34.94 
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32.13 
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30.00 
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28.001 


27.00 


26.50 


26.00 


25.50 


25.001 


24.50 


24.001 


23.50 


23.001 


22.50 


22.00 


21.501 


21.00 


20.501 


20.00 


19.50 


19.00 


18.501 


18.00 


17.501 


17.00 


16.50 


16.00 


15.501 


15.00 


14.50 


14.00 


13.50 


13.00 


12.501 
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11.00 
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11.00 
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11.00 
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11.00 
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6085.45 


6118.12 


6151.04 


6184.23 


6217.70 


6251.451 


6285.47 


6354.35 


6389.21 


6424.351 


6459.77 


6495.47 


6531.461 


6567.73 


6604.29 


6641.131 


6678.26 


6715.67 


6753.381 


6791.37 


6829.65 


6868.22 


6907.09 


6946.24 


6985.691 


7065.48 


7105.82 


7146.45 


7187.38 


7228.611 


7270.14 


7311.97 


7354.10 


7396.53 


7439.27 


7482.3ll 


7525.66 


7569.31 


7613.27 


7657.57 


7702.23 


7747.25 


7792.63 


7838.37 


7884.461 


7930.92 


7977.73 


8024.891 


8072.41 


8120.28 


8168.511 


27.29 


26.81 


26.56 


26.32 


26.07 


25.82 


25.57 


25.32 


25.07 


24.82 


24.57 


24.31 


24.061 


23.80 


23.55 


23.291 


23.03 


22.77 


22.52 


22.25 


21.99 


21.73 


21.471 


21.21 


20.94 


20.681 


20.41 


20.15 


19.88 


19.61 


19.35 


19.08 


18.81 


18.541 


18.27 


18.00 


17.731 


17.46 


17.19 


16.93 


16.67 


16.41 


16.16 


15.9 ll 


15.66 


15.41 


15.171 


14.93 


14.691 


1739985.13 


1739988.011 


1739989.23 


1739994.83 


1739995.79 


1739996.20 


1739997.18 


1739997.45 


1739997.68 


1739997.90 


1739998.09 


1739998.27 


1739998.56 


1739998.69 


1739999.09 


1739999.17 


1739999.24 


1739999.30 


1739999.46 


1739999.50 


1739999.54 
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11.00 


ll.OOl 


11.00 


11.00 


ll.OOl 


11.00 


ll.OOl 


11.00 


11.00 


11.00 


ll.OOl 


11.00 


11.00 


ll.OOl 


11.00 


11.00 


ll.OOl 


11.00 


11.00 


ll.OOl 


11.00 


11.00 


11.00 


ll.OOl 


11.00 


ll.OOl 


11.00 


11.00 


11.00 


11.00 


ll.OOl 


11.00 


11.00 


ll.OOl 
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ll.OOl 


11.00 


11.00 


ll.OOl 


11.00 


11.00 


11.00 


11.00 


ll.OOl 


11.00 


ll.OOl 


11.00 


11.00 


11.00 


11.00 


ll.OOl 


11.00 


11.00 
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8217.09 


8266.02 


8315.30 


8364.93 


8414.92 


8465.25 


8515.94 


8566.97 


8618.35 


8670.08 


8722.16 


8774.59 


8827.37 


8880.50 


8933.98 


8987.81 


9041.98 


9096.51 


9151.38 


9206.61 


9262.19 


9318.12 


9374.40 


9431.03 


9488.02 


9545.36 


9603.05 


9661.10 


9719.51 


9778.28 


9837.40 


9896.88 


9956.72 


10016.92 


10077.49 


10138.42 


10199.72 


10261.38 


10323.41 


10385.81 


10448.58 


10511.73 


10575.24 


10639.14 


10703.41 


10768.06 


10833.10 


10898.52 


10964.32 


11030.52 


11097.10 


11164.07 


11231.44 


14.45 


14.22 


13.99 


13.77 


13.54 


13.32 


13.10 


12.88 


12.67 


12.46 


12.25 


12.04 


11.84 


11.64 


11.44 


11.24 


11.05 


10.86 


10.67 


10.48 


10.29 


10.11 


9.93 


9.75 


9.58 


9.40 


9.23 


9.06 


8.89 


8.73 


8.56 


8.40 


8.24 


8.09 


7.93 


7.78 


7.63 


7.48 


7.33 


7.18 


7.04 


6.90 


6.76 


6.62 


6.48 


6.35 


6.22 


6.09 


5.96 


5.83 


5.70 


5.58 


5.46 


1739999.58 


1739999.61 


1739999.64 


1739999.67 


1739999.70 


1739999.72 


1739999.741 


1739999.761 


1739999.78 


1739999.79 


1739999.81 


1739999.82 


1739999.83 


1739999.851 


1739999.861 


1739999.88 


1739999.881 


1739999.90 


1739999.90 


1739999.91 


1739999.92 


1739999.92 


1739999.93 


1739999.93 


1739999.94 


1739999.94 


1739999.94 


1739999.95 


1739999.95 


1739999.95 


1739999.95 


1739999.96 


1739999.96 


1739999.96 


1739999.96 


1739999.97 


1739999.97 


1739999.97 


1739999.97 


1739999.97 


1739999.97 


1739999.97 


1739999.98 


1739999.98 


1739999.98 


1739999.98 


1739999.98 


1739999.98 


1739999.98 


1739999.98 


1739999.98 


0.000007 


0.000007 


0.000006 


0.000006 


0.000005 


0.000005 


0.000004 


0.000004 


0.000004 


0.000004 


0.000003 


0.000003 


0.000003 


0.000003 


0.000002 


0.000002 


0.000002 


0.000000 


0.000000 


0 . 000000 



0.000002 


0.000001 


0.000001 


0.000001 


0.000001 


0.000001 


0.000001 


0.000001 


0.000001 


0.000001 


0.000001 


0.000001 


0.000001 


0.000001 


0.000001 


0.000001 


0.000001 


0.000001 


0.000001 


0.000001 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 


o.ocoooo 


0.000000 


0.003000 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 
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11.00 


11.00 


11.00 


11.00 


11.001 


11.00 


11.00 


11.001 


11.00 


11.00 


11.001 


11.00 


11.00 


11.00 


11.00 


11.001 


11.00 


11.00 


11.00 


11.00 


ll.OOl 


11.00 


11.001 


11.00 


11.00 


11.00 


11.00 


ll.OOl 


11.00 


11.00 


11.00 


ll.OOl 
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Figure 3.0.5: Plot of Weight vs. Time 
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Figure 3.0.6: Plot of Inclination vs. Time 
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Figure 3.0.7: Plot of Total Velocity vs. Time 
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Figure 3.0.9: Plot of Thrust vs. Time 
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Figure 3.0.10: Plot of Thrust/Weight Ratio vs. Time 
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Figure 3.0.11: Plot of Y-Position vs. X-Position 
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4.0 Orbital Analysis 

The NASA ground rules for the Shuttle II specified two orbits into which the orbiter must be 
capable of launching: a 270 nautical mile orbit due east of Kennedy Space Center, and a 150 
nautical mile polar orbit To verify the feasibility of the proposed ASTS configuration, it 
was necessary to determine the types and magnitudes of the orbital maneuvers that would 
be necessary for orbit inserton. As a starting point, the output data for the final position and 
velocity from the launch simulation program was used. This data places the orbiter at an 
initial altitude of 78.1 nautical miles with a velocity of 25,962 feet per second. 


4.1 The Insertion Orbit 

The initial position and velocity of the orbiter place the orbiter originally at the perigee of 
an elliptical orbit. The semi major axis of this orbit, ai, was calculated using the Vis-Viva 

equation: . ^ ^ 

where ri is the sum of the altitude of the orbit and the radius of the earth (20925672.57 feet). 

The eccentricity of this first orbit, and the altitude and velocity at the apogee of this orbit 
were calculated using simple geometric relationships. The pertinent data on the first orbit 
is as follows: 


ei = .02473 
ai = 21,943,244 feet 
nap = 22,485,871 feet 
Vlap = 24,709 feet/second 


4.2 The 270 Nautical Mile Orbit 

In order to provide a minimum energy transfer between the insertion orbit and the 270 nmi 
orbit, an elliptical transfer orbit was used (see figure 4.2). 

The 270 nmi orbit has a radius of a2 = 22,567,272 ft. The semi major axis of the transfer 
ellipse, at, has a value of (ai +a 2)/2 = 22,255,258 ft. 


The velocities at periapsis and apoapsis of the transfer ellipse, as well as the velocity for 
the circular orbit were calculated using the Vis-Viva equation again. 


^ (t, ' a T ) 



7 







The following values were calculated for the velocities and the velocity changes needed for 
insertion into a 270 nmi orbit: 


Viap 24,709 ft/sec 

V t i 24,890 ft/sec 

dVi 181 ft/sec 


V t 2 24,799 ft/sec 

V 270 24,975 ft/sec 

dV2 176 ft/sec 

The total impulsice change in velocity for this orbital maneuver is 351 ft/sec. 


4.3 The 150 Nautical Mile Orbit 

The path of the elliptical orbit transfer orbit used for this maneuver is shown in figure 4.3. 

The velocity and radius of the 150 nmi orbit are found to be 

r 150 21837672 ft 

Viso ’ 25389 ft/sec 

Following the sane procedure as before in order to obtain the velocity changes necessary 
yields the following information: 

Vi a p 24709 ft/sec 

Vti 24678 ft/sec 

dVi -31 ft/sec 

Vt2 25419 ft/sec 

Vi 5 o 25389 ft/sec 

dV2 -30 ft/sec 


The total impulsive change in velocity for this orbital transfer is 
velocity changes would have to be made in a direction opposite to 
the shuttle. 


-61 ft/sec. These two 
the direction of travel 


of 
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Figure 4.3.1:150 Nautical Mile Polar Orbital Transfer 
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5.0 Stability and Control Analysis 

Stability and Control Analysis for the Advanced Space Transportation System was 
performed using geometric Shuttle data and empirical and graphical methods given by 
Etkin, Roskam and Smetana. The geometric parameters used to perform the analysis 
included physical data for the wing, the body, and the vertical tail. 

In order to determine stability and control derivatives, the geometric data for the Shuttle II 
must be used to determine aerodynamic parameters such as aspect ratio, wing taper ratio, 
and the mean aerodynamic chord. 

'A' is the aspect ratio and is defined to be equal to b^S; b is the wing span and S is the wing 
planform area. 'S' can be easily determined using the following formula: S=(c r +ct)(b/2), 
where c r is the wing root chord and ct is the wing tip chord. The ratio of the wing tip chord to 
the wing root chord is in itself and important physical parameter known simply as the 
taper ratio, X. The value of the mean aerodynamic chord (the wing chord which passes 
through the centroid of either half-wing) can be found by integrating the equation for the 
square of the wing chord: c 2 = {c r - [(c r -c t )/(b/2)]y) 2 over the span of half a wing and 
multiplying by 2/S. All of these aerodynamic parameters have been calculated and are 
found in the Executi ve Summary in Tables 3.2.1 and 3.3.1. ___ 

5.1 Determination of Clcc 

CLa is an extremely important stability derivative describing the change in lift of an 
aircraft due to variation in angle of attack. 

Roskam gives the following formula for determining CLa : 

CLa = KwbCLaw 

where K w b =1 - .25(d/b) 2 + 0.25(d/b) 

and CLaw = 2jtA/(2+[(A 2 B 2 /K 2 )(l+tan 2 Ac/2/B 2 )+4] 1/2 } 

In these equations: d = body diameter at root chord = 30 ft. 

b = wing span = 127.7 ft. 

A = wing aspect ratio = 2.8 

P = compressibility parameter = (1 - M 2 )^ 2 (for subsonic flight) 

K = ratio of actual wing section left curve slope to 2n 
Ac /2 = sweep of the wing half-chord line = 34° 


Since these equations are useful only for Mach numbers less than the critical Mach 
number, K is assumed to be approximately 1. Hence, the following values for CLa are 
obtained: 


M range P? 

0.0 - 0.3 1.0 

0.3 - 0.6 0.775 

0 .6-Mcr 0.415 


(/rad) 

2.97 

3.09 

3.33 


It is apparent that these values alone for the lift-curve slope (although correct) cannot be 
used to represent Cl<x for the ASTS. For this reason, methods of estimating CLa higher 
Mach numbers and varying Reynolds numbers (from Dynamics of Flight, by Etkin) were 
employed. 
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Etkin presents a graphical method by which the lift-curve slope of an aircraft may be 
obtained as a function of its geometrical parameters. Before incorporating Etkin s met , 
values for t/c and X (defined in Figure 5.1.1) were assumed to be 0.06 and 5 , respectively. 


Using Figure 5.1.1 the following section lift curve slopes were obtained: 


Bn 

lailflZlailoT 

iailQj/radl 

100 

0.87 

5.74 

10 7 

0.93 

6.14 

10 8 

0.95 

6.27 


It follows that these values for the actual lift curve slope are good approximations for 
transition both at the leading edge and at c/2. 

These values of (ai)o were then used in conjunction with Figure 5.1.2 to estimate the 
overall aircraft lift curve slope. Also, Etkin defines K as [(ai)om/2*]P and Ap as the 
compressibility sweep parameter equal to tan* 1 (tanAc/ 4 /p). 


Table 5.1.1. CLa of ASTS for M < Mcr 


M Range 

P 

Ab (°) 

Rn 

K 

pA/k 

pCLa/K </°) 

CLa(/rad) 

0.0 - 0.3 

1 

42 

io<J 

0.914 

3.06 

0.051 

2.67 

0.0 - 0.3 

1 

42 

i 0 r 

0.977 

2.81 

0.050 

2.80 

0.0-0.3 

1 

42 

lO 8 ^ 

0.998 

2.805 

0.049 

2.80 


0.88 

46 

10 s 

0.804 

3.06 

0.049 

2.57 


0.88 

46 

10 7 

0.860 

2.87 

0.048 

2.69 


0.88 

46 

10 8 

0.878 

2.81 

0.048 

2.74 

0.6 - Mcr 

0.644 

| 54 

10 * 

0.589 

3.06 

0.0451 

2.36 

0.6 - Mcr 

0.644 

54 

10 ' 

0.629 

2.87 

0.043 

2.41 

0.6 - Mcr 

0.644 

54 

10 8 

0.643 

2.80 

0.043 

2.46 


Note that these values of Cl« are slightly higher than those estimated using the method 
outlined by Roskam. 

Etkin presents a system for predicting Cl« for Mach numbers greater than the critical 
Mach n um ber using the graphical method shown in Figure 5.1.3 which is a function o 
M, X, and Ac/2- Results using this figure are listed in Table 5.1.2. 
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Table 5.1.2. CLa of ASTS for M > Mcr 


iMranae ! 

AtanAc /2 

A(M 2 -1)^ 

CLa (/rad) 

ekbiH 

1.89 

1.64 

4.14 

EEBEfli 

1.89 

2.97 

3.42 


1.89 

r 4.04 

2.72 i 


1.89 

5.02 

1 2.21 

tAAJvH 

L89] 

f 5.97 

P““ 1.82 

11 

1.89 

6.88 

1.62 


rrecucuon ui — - 0 - 

the lack of a proper estimation technique. 

5.2 Determination of the Neutral Point, hn 

The overall aircraft aeutral point, hn, i. defined a. the point where the pitching moment it 
invariant with angle of attack. It i» further defined geometncally by Figure 5.2.1. 

Roskam gives the following relationship for determining the neutral point . 

hn = hnw + Ahnfl (for tail-less configurations) 

where hnw = ki[(XacVc r )-k2] 

and AhnB - - (dM/da)/(qScCLaw) 

where dM/da = (q/36.5) S Wf 2 (Xi)de/daAXi 

Note that Roskam defines both the neutral point, hn, and the aerodynamic center, Xac/c, as 
the same quantity. 

Using the relationships for Ahn B , it was found that the body had little or no effect on the 
overall aircraft neutral point Therefore, determination of the neutral point is reduced to 
the following form: 

hn = ki[(XacVc r )-k 2 ] 

Xac'/cr ki and k 2 are found using Figures 5.2.2, 5.2.3, and 5.2.4, respectively Using 
these figures together with known geometric data and the relationship for hn, the following 
results were obtained. 


Table 52.1. Neutral Point of the ASTS for M < Mcr 


M range 

P 


bi 

k 2 

hn 

0.0 - 0.3 

1.0 


1.43 

0.39 

0.27 

0.3 - 0.6 


0.59 

L43 

0.39 

0.29 

0.6 - Mcr 

0.64 

0.60l 

1.43 

0.39 

0.30 
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These values for hn show that as the Mach is increased, the neutral point begins to move 
gradually aft. This finding agrees with most aerodynamic data on similar aircratt. 

It is important to determine the characteristics of a flight vehicle in as many of its flight 
regimes as possible. Therefore, methods given by Etkin were used to find values of hn for 
Mach numbers greater than the critical Mach number, 

Etkin gives a graphical method shown in Figure 5.2.5 for determining hn as a function of 
A, MA and A d2- Using these variables and Figure 5.2.5, the following data was obtained. 


Table 5A2. Neutral Point of the ASTSfbrO< M<2£ 


M range 


■HBSEIS&Syi 

hn 

0.0 - 0.3 



0.28 


1.9 

2.5 

0.29 

lo.6 - Mcr 

1.9 

1.8 



1.9 

1.6 


11.3 -1.6 

r i-9 

30 



1.9 

4.0 

o.4a 

HUM 

1-9 

5.0l 

0.48 

E OE23H 

1.9 

6.0 

r ~ 0.49 


1.9 

S 6.9 

0.49 


Again, these numbers show that (as expected) the neutral point starts out relatively close to 
the quarter-chord for low subsonic Mach numbers and moves aft approaching the halt- 
chord at higher supersonic Mach numbers. 


5.3 Determination of Cmo 

C mo is an extremely important stability derivative which dictates whether or not a 
configuration has positive or negative static stability- Cmo is determined through the 
following equation: 


Cmo = CLct (b - hn) 


where h represents the non-dimensional location of the aircraft center of gravity. For ® 
configuration to be statically stable, Cmo must be negative; this requires that the quantity 
(h - hn) must be negative. Xcg was determined by considering the location and weight of 
the components of the ASTS. The distance between the center of gravity and the leading 
edge of the mean aerodynamic chord (non-dimensionalized by dividing by the mean 
aerodynamic chord) is defined as h. Upon completion of these calculations, it is found that 

h = 0.136. 


Hence, Cmo = Cj/x (0.136 - hn) 


5.4 Cia of Vertical Tail 
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The lift curve slope of the vertical tail is determined by the geometric characteristics of the 
vertical tail and is influenced by the change in flow characteristics due to the presence of 
the wing and body combination. Figure 5.4.1 shows the relationship between (ap)B/ai for a 
vertical tail on a body of circular cross section. Using a value of 1 for D/h, (aF)B/ai is 
approximately 2.9. 

Values ofcEafor the vertical tail are determined using the same procedure as before for the 
wing. When-this value is obtained it is multiplied by the factor (aF)B/ai to correct for its 
locationmthe flow field relative the the wing and body. It is assumed that the vertical tail 
has approximately the same values of t/c and t as the wing so that both the wing and 
vertical tail have the same values for (ai) 0 . Using the geometric parameters for the 
vertical tail and Figures 5.1.2 and 5.1.3 the following is found. 


Table 5.4.1. Corrected Clo values for the ASTS Vertical Tail 


M range 

Bn 

Cta 

Cl a corrected (rad* 1 ) 

0.0 - 0.3 

10 6 

2.93 

8.20 

0.0 - 0.3 

10 7 

3.08 

8.60 


10** 

3.14 

8.80 

lo.3-0.6 

”” —W 

2.62 

7.30 

EEBiXifli 

id 1 

2.74 

7.70 


10^ 

2.75 

| 7.70 

0.6 - Mcr 

10** 

2.36 

6.60 

0.6 - Mcr 

id 1 

P " 2.46 

6.90 

0.6 - Mcr 

- W 

2.46 

— — 6 9Q 


— 

' 4.58 

” 12.80 

mmsm 

... 

4.35 

12.201 


Determination of Cut corrected for the vertical tail beyond a Mach number of 
approximately 1.6 was not possible. 


5.5 Determination of C| P 

Cio is the change in rolling moment due to a variation in rolling velocity. Ci p is often 
referred to as the rolling moment damping derivative. Etkin provides a graphical method 
by which Ci p can be calculated based on the following parameters: P, KA*,and Ae (which 
is defined as tan* 1 [(l/pltanAc/^). Etkin’s method is presented here in Figure 5.5.1. Table 
5.5.1 contains the values of Ci p found using Etkin’s method. 
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Table 5.5.1. Cip values ofASTS for M < Mcr 



C]p does not vary greatly with Mach number, and is always negative. 

5.6 Determination of Cip 

Cm represents the change in rolling moment due to a change in sideslip angle, p. Qp is 
primarily governed by the wing dihedral angle, the body and the vertical tail For this 
reason, Cip is often referred to as the effective dihedral derivative. Etkin present a 
graphical method of determining Cip based on the effects of the wing, the body and the 
vertical tail. This method utilizes Figure 5.6.1. Using this figure, it was found that 
Cipw = *.32 Cl. 

The effect of the body on C]p is found using the equation: 

(ACip)B = 1-2 (A) 1/2 Zw (h+wVb 2 

where Zw = vertical distance of c/4 point below fuselage centerline = 11.0 ft. 
h = average fuselage height at wing root = 28.0 ft. 
w = average fuselage width at wing root = 30.0 ft. 

Using this information: (ACip)B = 0.08. 

The effect of the vertical tail on Cip is found using the equation. 


(ACip)p = -ap (1 - do/dp)(Vp/V)2 Sp Zf/(S b) 

where Zp = vertical distance between vertical tail AC and aircraft CG = 41.0 ft. 
Sp = Vertical Tail planform area = 603 ft 2 . 

Using this equation under the assumption that da/dp is negligible and V F /V is 
approximately equal to 1: (ACip)F = -0.033 ap 


Combining components: Cip= -0.32 Cl + 0.08 - 0.033ap 

This equation holds for .til Mach numbers below the critical Mach number. 
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5.7 Determination of Cir 


C lr represents the change in rolling moment due to yawing velocity, or simply rolling 
moment due to yaw. Etkin divides Ci r into two components: Ci r due to the wing an ir 
to the vertical tail. Using Figure 5.7.1, (C] r ) B = 0-23 C L• 


due 


The contribution of the vertical tail on Ci r is shown through the following formula. 


(Cjr )f = aF Sf Zf [z +(lf/b)]/(S b) 

where If = distance between vertical tail AC and the aircraft CG = 48.0 ft. 


Hence, (Q r )F = 0-079 aF 


Combining components: C] r = 0.23 Cl + 0.079 ap 

This equation holds for all Mach numbers below the critical Mach number. 


5.8 Determination of C n p 


C n B is defined to represent the change in yawing moment due to sideslip angle, p. C n p is a 
very important stability derivative since it alone determines whether or not an aircraft has 
positive weathercock stability; for positive weathercock stability, C n p must be positive. 


Figure 5.8.1 shows the variation of C n (j/C L 2 as a function of wing aspect ratio, A, and wing 
quarter-chord sweep angle, A. 


Using Figure 5.8.1, C n p = 0.07 Cl 2 

5.9 Determination of C n p 

C np is a cross derivative defined as the change in yawing moment due to rolling velocity. 
Etkin presents the following equation for determining Cn p : 

C np = [(AC np )l/ClJ Cl + KACnpWDOla] (Cdo)<x 

where (AC np )l/CL and (AC np )2/(CD0)a are found from Figure 5.9.1 

From Figure 5.9.1, (AC np )l/CL = -0.17 and (AC^Cdo* = 19 > *> that C np may be written 

3 S 

Cnp = -0.17 Cl + 19 (Cdo)cl 

5.10 Determination of Cnr 

Cnr is defined to be the change in yawing moment due to yawing velocity and is often 
referred to as the damping in yaw derivative. Etkin presents the following equation for 
determining C n ri 

Cnr = KAC nr )i/CL 2 ]CL 2 + [(AC nr )2/CDo]CDO 
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where [(ACnrh/CL 2 ] and [(AC nr )2/CD0] are found from Figure 5.10.1 

Figure 5.10.1 reveals that [(ACnrh/CL 2 ] = -0-02 and [(ACnr^CDo] = -0.48, so that C nr may be 
written as 

Cm. = -0.02 Cl 2 - 0.48 Cdo 

5.11 Determination of CMq 

C mq is defined as the variation in pitching moment due to pitch velocity or simply as the 
damping in pitch. Smetana gives the following relationship for Cmq : 

C mq * - 2x* lx'I Clo/c 2 

where x' = distance from CG to wing quarter-chord (positive for forward CG) = -6.97 ft. 

After substituting the proper geometric and stability parameters, C m q can be written as 
Cmq = 0.0372 CLa 

5.12 Determination of Ciq 

CLq is defined as the change in lift coefficient due to pitch velocity. Smetana gives the 
following relationship for CLq • 

Cj^q — 2 x Cjja/c 

Upon substitution of known geometric parameters, CLq is written as 
CLq = -0.2727 Clo 

5.13 Determination of Cd<x 

CDa is defined as the change in drag coefficient due to a change in the angle of attack. 
Smetana gives the following relationship for CDa: 

CDa = (dCDO/doO + [2 Ci/frA)] Clo 

Since (dCoo/da) is usually small, CDa may be written simply as [2 Cl/(*A)] CLa- 

Therefore, CDa may be expressed as 
Cdo = 0.2274 Cl CLa- 

5.14 Determination of Cyp 

Cyft is defined as the change in side force,y, due to a change in sideslip angle, P- Smetana 
breaks Cyp up into 3 components: that of the wing, body, and vertical tail. 

(Cyp)w = Cl 2 (6 tan A sin A)/[MA + 4 cosA)] 

Upon substitution, (Cyp) w = 0.2232 Cl 2 
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(C y p)B = -ki (Clo)b (Sb/S) 

where Sb = body reference area = 3733 ft^ 

ki = interference factor determined from Figure 5.14.1 


Figure 5.14.1 gives ki as a function of Z*/(d/2) where Zw is the distance from 1body 
centerline to quarter-chord point of exposed wing root chord (positive for quarter-chord 
point below body centerline) and d is the max body height at wing/body intersection. 


Zw=ll-0ft. 
d = 25.0 ft 


From Figure 5.14.1, ki = 1.44 

Smetana gives the following values of (Clo)B ; 

(Cl<i) circular body = 0.0525 rad' 1 
(Cl<x) rectangular body = 0.1253 rad' 1 

Since the ASTS configuration is a combination of these two shapes, let (Cl<x)b - 0.08 rad 

Now that all values needed are known, (C y p)fi = -0.074 

(Cyp)F = -k (Clo)f (1 + da/d[3) Sp/S 

k is found from Figure 5.14.2 and is a function of (W(2ri)]. 

by = vertical tail span = 29.0 ft. 

ri = fuselage diameter in tail area = 30 ft. 

From Figure 5.14.2, k = 0.75 

If it is assumed that do/dB is negligible then (Cyp)F may be expressed as 
(Cyp)F = -0.078 aF 

Combining all components; Cyp = 0.2232 Cl^ - 0.074 - 0.078 aF 


5.15 Determination of C yp 

Cyp represents the change in side force due to rolling velocity. Smetana expresses C yp as 
the sum of the contributions due to the wing and the vertical tail. 

(Cyp ) w = Cl ([(A + cosA) tanMA + 4 cosA)] + 1/A) 

Substitution into (Cyp )w yields (Cyp )w = 0.91 Cl 

(Cyp )f = -2 aF Sf Z F /(S b) 

All of these parameters have been defined previously so that (Cyp )F = -0.067 aF 
Combining these two; Cyp = 0.91 Cl - 0.067 aF 
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5.16 Determination of Cyr 

Cyr is defined as the change in side force due to yawing velocity. Smetana expresses Cyr 
asthe sum of the contributions due to the wing and the vertical tail. 


(Cyr)w = 0.143 Cl * 0.05 


(Cyr)p = -2 If (Cyr)F/b = -0.752 (Cyr)F 
So that, Cyr = 0.143 Cl - 0.05 - 0.752 (Cyr>F 


5.17 Dynamic Stability Matrices 

Once the proper stability derivatives have been determined, it is possible to use them to 
investigate the dynamic stability of the configuration. Presented here are the longitudinal 
and lateral stability matrices corresponding to the non-dimensional, controls fixed, 
linearized equations of motion for a rigid aircraft. 


Longitudin al Eouations: 

Dx = Ax where x = (u a q 0) T 

A(l,l)»(Cxu + 2 Clo taneoV2n 
A(l,2) = Cxo/2p. 

A(l,3) = 0.0 

A(l,4)»-CL0/2p 

A(2,l) * (Czu * 2 CloV(2p - Czo) 

A( 2 , 2 ) = Czo/(2|X - Czo) 

A(2,3) = (2p + CzqV(2p - Czo) 

A(2,4) = -Clo tan©o/(2p - Cz a) 

A(3,l) = (C M u + Cmo A(2,l))/iB 
A(3,2) = (CMa + C mo A(2,2)ViB 
A(3,3) = (Cftfq + Cmo A(2,3))/iB 
A(3,4) = C ma AC2,4ViB 
A(4,l) = 0.0 
A(4,2) = 0.0 
A(4^) = 1.0 
A(4,4) = 0.0 


where p = m/(p S c/2) 

iB = iyy/tp s (c/2) 3 ] 

Cxa = Clo * CDa 
Cza = *Clo * Cdo 
Czu = * (M 2 /l-M 2 ) C lo 
CZq = -CLq 
CDm = CDu/M 

Cxu = -2 (Cdo + Clo tanQo) - M CDm 
Czo = 0.0 
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T.lateral Equations; 

Dy = By where y = (B p r <J>)^ 

B(l,l) = Cyp/2p 
B(l,2) = Cyp/2n 
B(l,3) = -(2 \i - CyrV2|i 
B(l,4) = Clo^IA 
B(2,l) = Gc Cip + iE C n pVA 
B(2,2) = (ic Ci p + i£ C np )/A 
B(2,3) = Gc Clr + iE CnrVA 
B(2,4) = 0.0 

B(3,l) = Ga Cnp + iE QpVA 

B(3,2) = Ga C np + iE Cl p )/A 

B(3,3) = GACnr + iEClryA 

B(3,4) = 0.0 

B(4,l) = 0.0 

B(4,2) = 1.0 

B(4,3) = tan6 0 

B(4,4) = 0.0 

where p = m/(p S h/2) 

A = iAiC*iE 2 
iA = Ixx/[p S (b/2) 3 ] 
ic = Izz/Ip S GV2) 3 ] 
ig = Ixz/[p S (h/2) 3 ] 

Once the matrices A and B are found, their Eigen Values can be determined. These Eigen 
values reflect whether or not the system is stable. For a system to asymptotically stable, all 
of its Eigen values must have negative real parts. 

The matrices A and B were calculated using a computer program written in FORTRAN. 
The program and results of the investigation of the dynamic stability for landing 
(assuming 0 O = 10° and U 0 = 300 ft/sec and L = W) are listed on pages 66 and 71. 


5.18 Determination of C|§ a 

Ci8a describes the variation in rolling moment with change in aileron deflection; often, 
this s tab ility derivative is referred to as the aileron power. Smetana gives the following 
equation for Cis®: 

Ab 

Ci 5 a = [2CLaX/(Sb)]\ cydy 

where c is the equation for the wing chord as a function of wing span, y. 


c = c r - [(c r -c t )/(h/2)] y = 73.125 - 0.863 y 

and a and b represent the spanwise distance for the body centerline to the inner and outer 
aileron points, respectively. 



a = 27.0 ft b = 58.5 ft 

x is a correction factor which is a function of the aileron chord to wing chord ratio, t is 
determined from Figure 5.18.1. Since the ratio of aileron chord to wing chord is not 
constant over the span of the aileron, the average value was used to determine x. 

(c a /c w )avg = 0.33 t=0.57 

Upon integration and substitution of the proper values, 

C] 5ct = 0.071 CLa. 

5.19 Determination of Cisc 

Cj 5 r is the variation in rolling moment coefficient with change in rudder deflection. 
Smetana gives the following formula for estimating C]gr • 

Qfir = ap t Sp Zp/(S b) 

Here, t is a function of the rudder area to vertical tail area ratio. 

Sr = 157.5/f0) Sf = 603 ft2 

Hence, S^F = 0.261 so that from Figure 5.19.1, x = 0.48. 

Substitution into the equation for C]gr results in 
CiSr = 0.0158 ap. 

5.20 Determination of C n sr 

Cn5r is the variation in yawing moment coefficient with a change in rudder deflection and 
is often referred to as the rudder power. Smetana gives the following relationship for C n 5r 
C n Sr = -*F * Sf lf/(S b) 

As with CiSr, the value of x is determined from Figure 5.19.1 and has a value of 0.48 so that 
Cn5r may be expressed as 

Cnfir = -0.0187 ap 
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AE448 CG Locations/1 



A 

B 

C 

D 

1 

Component 

X-Location 

Waiaht 

Y-Moment 

2 

Winer 

116.60 

28.300 

3,299,780 

3 

Vertical Tail 

140.00 

3,200 

448,000 

4 

Forward Body 

28.90 

27,000 

780,300 

5 

Mid Body 

80.00 

23,600 

1,888,000 

6 

Aft Body 

128.00 

13,200 

1,689,600 

7 

Main Enaines 

137.00 


3,123,600 

8 

OMS 

138.00 


r 510,600 

9 

RCS 

138.00 

1,600 

220,800 

m 

Control Systems 

48.00 

4,000 

192,000 

Vfl 

Avionics 

22.00 

3,100 

68,200 

m 

Personnel 

35.00 

1,600 

56,000 

m 

Stores 

1 42.00 

600 

25,200 

m 

Encrine Accessories 

125.00 

6,200 

775,000 

m 

Hvdroaen Tanks 

64.00 

24,000 

1,536,000 

m 

Oxvaen Tank 

121.00 

11,000 

1,331,000 

eh 

TPS 

96.30 

6,100 

587,430 


Totals 


180,000 

16,531,510 

m 

Centroids 

91.84 
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a K> CM 


V 

V. 



f' He'!<;'< CONTROI S FIXED LONGITUDINAL.. AMD LAIERhL 
C *** DYNAMIC STABILITY AND CONTROL PROGRAM 
C *#* FOR THE SECOND GENERATION SPACE 
C #** TRANSPORTATION SYSTEM. 

C ### WRITTEN BYi PHILIP BENEFIELD & DEL JOHNSON 


C 


C 


C 


DOUBLE 

DOUBLE 

DOUBLE 

DOUBLE 

DOUBLE 

DOUBLE 

DOUBLE 

DOUBLE 

DOUBLE 


PRECISION G.S » W»US.CBAR .UO .RHO .THETAO 
PRECISION B . A » C . E . C L 0 . CDO . CMU . CM A 
PRECISION CMAD»CMQ.CLA»CLQ.CDA.CDU 
PRECISION CL.B . CLP v CLR . CNB» CNP.CNR 
PRECISION CYB .CYP .CYR .LLONG.MULGNG.IBLONG 
PREC I SION CXA . CZA . SOUND » M.CZl.1.CZQ . CDM 
PRECISION CXU » CZAD. ALONG(4.4) »LLAT 
PRECISION IALAT .ICLAT»IELAT»MULAT 
PRECISION DEL.ALAT(4*4) 


OPEN(UNIT=9»FILE='5TAB1IN.DAT'»ST ATUS= ' OLD') 

• OPEN < UNIT*10.FILE*'STAB1OUT.DAT".STATUS*'UNKNOWN') 


READ (9. & ) G.S.W.US. CBAR . UO » RHO » THET AO 

READ(9»#) B.A.C.E.CLO .CDO 

READ(9»#) CMU.CMA.CMAD.CMQ 

READ(9. *) CLA .CLQ .CDA»CDU 

"READ ( 9»*) CLB.CLP . CLR 

READ(9»*) CNB »CNP . CNR . CYB .CYP .CYR 


00 

0 

* 

c 



FORMAT (IX»A40 .' = / »F16.6) 

FORMAT(////////.1OX t 'GEOMETRICAL * INERT IAL 
' CHARACTERISTICS: ',//) 


X EQUILIBRIUM 


y 


WRITE <10.20) 
WRITE( 10» 300) 
WRITE( 10. 300) 
WRITE( 10 . 300) 
WRITE<10.300) 
WRITE(10.300) 
WRITE <10. 300) 
WRITE(10.300) 
WRITE(10.300) 

WRITE <10.300) 
WRITE(10.300) 
WRITE(10»300) 
WRITE(10>300) 
WRITE(10. 300) 
WRITE(10.300) 


'GRAVITATIONAL ACCELERATION (FT/SEC 
'WING AREA (FT^Z)'.S 
'WEIGHT (LB)'.U 
'WING SPAN (FT)'.WS 

'MEAN AERODYNAMIC CHORD (FT)'.CBAR 

'AIRSPEED (FT/SEC) ' . UO 

'AIR DENSITY (SLUGS/FT”3)'»RHO 


INITIAL THETA (RAD) 


' IYY 

(SLUG-FT-2)'.B 

' IXX 

(SLUG-FT~2)'.A 

' IZZ 

(SLUG-FT"2)'.C 

' IXZ 

(SLUG-FT-2)'.E 

'CLO 

(-) ' .CLO 

'CDO 

(-) ' . CDO 


THET AO 


2) ' . G 


FORMAT ( //// :> 1 OX > ' NON-DIMENSIONAL ST ABILITY 
WRITE(10.30) 


DERIVATIVESt ' » 


WRITE(10.300) 
WRITE(10.300) 
WRITE(10.300) 
WRITE(10.300) 
WRITE(10.300) 
WRITE(10.300) 
WRITE(10.300) 
WRITE (10.300) 


CMU (-) ' » CMU 


CMA (-> ' » CMA 

Q?X 

CMAD (-)'» CMAD 

OF 1 

CMC (-) ' » CMQ 

CLA (-) ' . CLA 

* 

CLQ (-) ' . CLQ 

CDA (-) ' . CDA 

CDU '.CDU 



a nn 
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FOR MAT <//// p10X?'LATERAL 


DI RECTI ON At. 


WRITE(10?40) 
WRITE(10 ? 300) 
WRITE (10 ? 300 ) 
WRITEC10*300) 
WRITE(10 ? 300) 
WRITE(10 ? 300) 
WRITE(10 ? 300) 
WRITE < 10 ? 300) 
WRITE(10 ? 300) 
WRITE(10 v 300) 


' CLB (/RAD)'»CLB 
'CLP </RAD)'?CLP 
' CLR </RAD)'?CLR 
'CNB </RAD)'?CNB 
'CNF (/RAD)' f CNP 
'CNR </RAD)'rCNR 
'CYB </RAD)'?CYB 
'CYP (/RAD)'?CYP 
'CYR (/RAD)'?CYR 


derivatives: 


/) 


5(c^C5{t CALCULATIONS FOR LONGITUDINAL MOTION 


C 


LL0NG=CBAR/2♦0 

MULONG=<W/G)/(RHOTSTLLONG) 

IBL0NG=B/(RH0#S*(LL0NG#*3) ) 

CXA=CLO-CDA 
CZA=-CLA-CDO 
S0UND=1116*4 
M=UO/SOUND 

CZU=-((M*M)/<1.0-M*M)>*CLO 
CZQ=-CLQ 

CXU=“<2^0T <CDO+CLOTTAN(THETAO) ) )-<M*CDM) 
CZAD=0.0 


:*LONG (1 ? 1 > = 
ALONG(1»2) = 
ALONG <1 ? 3) 
ALONG(1 ? 4) 
ALONG( 2 ? 1 ) 
ALONG < 2 p 2 ) = 
ALONG < 2 ? 3) 
ALONG (2 p 4) 
ALONG < 3 p1) 
ALONG(3 p 2) 
ALONG( 3 p 3) 
ALONG(3 ? 4) 
ALONG<4?1) 
ALONG(4?2) 
ALONG(4v3) 
ALiJN(j v 4 ') 


(CXU + 2.OTCLOTT AN< THETAO))/(2.0*MUL0NG) 
CXA/(2.0*MUL0NG) 

0.0 

-CL0/(2.0*MUL0NG) 

(CZU-2.OTCLO)/(2.OTMULONG-CZAD) 
CZA/(2.0*MUL0NG-CZAD) 

(2♦OTMULONG+CZG)/< 2♦OTMULONG-CZAD) 

■ (-CLO*TAN(THETAO))/< 2.0*MUL0NG-CZAD) 
~(CMUTCMAD*AL0NG(2p1))/IBLONG 
=(CMAtCMADTALGNG<2p 2))/IBLONG 
ss < CMQFCMADTALONG ( 2 p 3) ) /IBLONG 
=CMADTALONG(2 ?4)/IBLONG 
= 0 . 0 
= 0.0 
= 1.0 
'= 0.0 


c *:!<* CALCULATIONS FOR LATERAL MOTION 

C 

LLAT=WS/2.0 

MULAT=<W/G)/(RKOTSTLLAT) 

IALAT-A/ ( RH0T.FaLLATr.K3 ) 

ICLAT=C/(RHCTSTLLATTT3) 

I.ELAT=E/ ( RH0TSTLLATTT3 ) 

DEL=<IALAT#ICLAT)-IELAT#*2 

C 



c 


10 

c 


15 

C 

100 


101 


200 

201 

c 


ALAT < 1 y 1) ~CYB/ ( 2. OY-MULAT ) 

AL AT < 1 9 2 ) ~CYP/ ( 2 * OfcMULAT ) 

ALAT (1 3) «- ( 2 * OYMULAT-CYE>/(?-♦ 0#MULAT ) 
ALAT (1» 4) -CL.O/ < 2 * 0*MULAT) 

ALAT(291>=(ICLAT*CLB+IELAT*CNB>/DEL 
ALAT(2 9 2) = (ICLAT*CLP+IELAT*CNP)/DEL 
ALAT ( 2 9 3) - ( XCLAT*CLR +1EI..AT*CNR) /DEL 
ALAT(294)=0♦0 

ALAT(3 9 1) = (IALAT#CNB+IELAT*CLB)/DEL 

ALAT(3 9 2) = (IALAT*CNP+IELAT*CLP)/DEL 

ALAT(3 9 3) = (IALAT#CNR+IELAT*CLR)/DEL 

ALAT(39 4)=0*0 

ALAT(4 y1)=0♦0 

ALAT(4 9 2)=1♦0 

ALAT (4 9 3 )=TAN(THETAO) 

- ALAT ( 474F-C/ 7tT~ 


WRITE(10 9 100) 

DO 10 J1 9 4 

WRITE( 109200) ( ALONG ( J 9 K) •> K=1.94) 

IF(J♦NE♦4) WRITE (109201) 

CONTINUE 


WRITE(10 9 101) 

DO 15 J=l94 

WRITE( 109200) (ALAT(J 9 K)9 K = l94) 

IF ( J ♦NE. 4) WRITE (109201) 
CONTINUE 


FORMAT(//////9 25X 9 ' THE STABILITY MATRIX CORRESPONDING TO ' 
THE NONDIMENSIONAL 9 CONTROLS FIXED 9 LINEARIZED t ' */* 
LONGITUDINAL EQUATIONS OF MOTION FOR' 9/9 

'A RIGID AIRPLANE' 9 //) , ... , 

■<////»25X»'THE STABILITY MATRIX CORRESPONDING TO 9 / 
jcrrnuAi . r.nNTROLS FIXED 9 LINEARIZED 9 ' 9 / 9 


9/9 


20X9 
25X9 
35X 9 ' 
FORMAT 
20X9 
23X9 
35X 9 
FORMAT 
FORMAT 


THE N 0N DIMENSION AL 9 CONTROLS FIXED 9 


LATERAL 
A RIGID 
(5X9'I 

(3X9 


EQUATIONS 
AIRPLANE' 9 //) 

'9 4(F12 ♦ o 9 6X )9 ' 
974X9'i') 


OF MOTION FOR' 9/9 


CLOSE(9) 
CLOSE(10) 
STOP 
END 


\ 
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GEOMETRICAL * INERTIAL 


& EQUILIBRIUM 


CHARACTERISTICS 


a ♦ 


GRAVITATIONAL ACCELERATION (FT/SEC 2) 

WING AREA (FT~2) - 
WEIGHT (LB) = 
WING SPAN (FT) = 
MEAN AERODYNAMIC CHORD (FT) = 
AIRSPEED (FT/SEC) = 
AIR DENSITY (5LUGS/FT~3) - 
INITIAL THETA (RAD) * 
IYY (SLUG-FT~2) = 
IXX (SLUG-FT"2) = 
IZZ (SLUG-FT"2) = 
IXZ (SLUG-FT"2) = 

CLO (-) = 

CDO (-) = 


32♦174000 
5818* 000000 
180000,000000 
127.692000 
51.120000 
300.000000 
0.002377 
0.174500 
12490393.000000 
1686164,000000 
13399266.000000 
795369.000000 
0.418000 
0.050000 


non-dimensional stability derivatives: 

CMU (-) = 

CMA (—) = 

CM AD (-) = 

CMQ (-) = 

CLA (-) = 

CLQ <-) = 

CDA (-) = 

CDU (-) = 


0.000000 

-0.403200 

0.000000 

0.104200 

2.800000 

-0.763500 

0.266000 

0.000000 


lateral-directional derivatives: 


CLB 

(/RAD) 

- 

CLP 

(/RAD) 

* 

CLR 

(/RAD) 

- 

CNB 

(/RAD) 


CNP 

(/RAD) 

= 

CNR 

(/RAD) 


CYB 

(/RAD) 


CYP 

(/RAD) 

- 

CYR 

(/RAD) 

- 


—0.338000 
-0.225000 
0.775500 
0.012000 
-0.070000 
-0.030000 
— 0 ♦ 706000 
-0,195800 
0.514200 




-0*00315899 
-0 * 02743687 
0*00000000 
0.00000000 


-0.05570903 

-0.74066202 

-0*04074188 

0.00000000 


TH 


n c- 


THE STABILITY MATRIX CORRESPONDING TO 
NONDIMENSIONALs> CONTROLS FIXED * LINEARIZED s- 
LONGITUDINAL EQUATIONS OF MOTION FOR 
A RIGID AIRPLANE 


0.00480166 

-0.09003117 

-0.00745438 

0.00000000 


0,00000000 

1.02411888 

0.00192645 

1.00000000 


- 0 * 0132045 / 
-0 . 00232787 
0 , 00000000 
0,00000000 


THE STABILITY MATRIX CORRESPONDING TO ^ 

THE NONDIMENSIONAL i CONTROLS FIXED» LINEAR!ZED» 
LATERAL EQUATIONS OF MOTION FOR 
A RIGID AIRPLANE 


-0.01545018 

-0.50320992 

-0.04867205 

1.00000000 


-0.95942552 

1,69903754 

0.09279544 

0.17629303 


0.03298353 

0.00000000 

0.00000000 

0.00000000 
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Figure 5.1.1. Lift curve slopes for two-dimensional incompressible flow. 






, per degree 



Figure 5.1.2.Lift-curve slope for swept and tapered wings 
number. 


at speeds below critical Mach 



per radian 



Figure 5.1.3. Lift-curve slope for swept and tapered wings at supersonic speeds. 
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Figure 5.2.2. Aerodynamic 



snter Location of Lifting Surfaces. 







0 .2 A . A *8 l-O 


Figure 5.2.3. Aerodynamic Center Transformation Constant Ki 



w 

w 








































Figure 5.2.5. Chordwise position of the mean aerodynamic center of swept and tapered 
wings at high speeds expressed as a fraction of the mean aerodynamic 


chord. 



<«F>» <*y>»r 



Figure 5.4.1. Lift-curve slope for single fin and rudder on a 


body of circular cross section 





A t , degrees 

Figure 5.5.1. C] p for straight-tapered wings. 
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*0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 


Figure 5.6.1. Cip for straight-tapered wings with no dihedral. 








Figure 5.7.1. Charts for estimating C,, for subsonic incompressible flow. 


A. degrees A. degrees 













Figure 5.9.1. Variation of (AC n ph/CL 

incompressible flow. 


and (AC np )2/(CD0)a with aspect ratio for subsonic 
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Figure 5.10.1. Charts for estimating (ACnrVCL 2 and (ACnr)2/C D 0 


a C. 













Figure 5.14.2. Values for k as a function of the ratio of vertical tail span to fuselage 
diamter in the tail region. 




Figure 5.18.1.Values for x as a function of aileron chord to wing chord ratio. 



Figure 5.19.1. Values for x hs a 


function of rudder area to vertical tail area ratio. 
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